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Abstract 

An approximate hadronic symmetry based on spin and flavor in- 
dependence and broken by spin and mass dependent terms is shown 
to follow from QCD. This symmetry justifies the SU(6) classification 
scheme, but is more general in allowing its supersymmetric extension 
based on a diquark-antiquark symmetry. It will be shown that the 
same supersymmetry is also implied in the skyrmion type effective 
Lagrangian which could be extracted from QCD. Predictions of the 
Skyrme model is improved by using different realizations of the chiral 
group. 



PACS numbers: 12.40.Aa, 12.40.Qq, 11.30.Pb 



f Work supported in part by DOE contracts No. DE-AC-0276 ER 03074 
and 03075, and PSC-CUNY Research Awards. 



1 



Introduction 



Inside rotationally excited baryons, QCD leads to the formation of diquarks 
well separated from the remaining quark. At this separation the scalar, 
spin independent, confining part of the effective QCD potential is dominant. 
Since QCD forces are also flavor independent, the force between the quark 
(q) and the diquark (qq) inside an excited baryon is essentially the same 
as the one between q and the antiquark q inside an excited meson. Thus, 
the approximate spin-flavor independence of hadronic physics expressed by 
SU(6) symmetry gets extended to the SU (6/21) supersymmetry through 
a symmetry between q and (qq), resulting into the parallelism of mesonic 
and baryonic Regge trajectories. Various aspects and implications of this 
approximate effective supersymmetry and its breaking are discussed in the 
first section. 

In the next section we discuss the Skyrmion approach to effective hadronic 

Lagrangians. In the early 1960 's Skyrme^ proposed a model in which 

baryons can be understood as solitons in an effective Lagrangian for bosons. 

[of 

He added a fourth term to the non-linear sigma models and proposed that 
the topological soliton solution of the meson field could be identified as a nu- 
cleon, and that the winding number is the baryon number. This pioneering 
work of Skyrme was largely ignored until two decades later when Pak and 
Tze^, and Gipson and Tze^ re-expressed Skyrme's ideas in the modern 
context of QCD and studied their implications for weak interactions. In an 
attempt to verify one of the Skyrme's conjectures that the winding numbers 
are the baryon numbers, Balachandran and his collaborators^ considered a 
system of interacting quarks and solitons. Using some results of Goldstone 
and Wilczek©, they showed that, in the presence of the classical soliton 
field, the Dirac sea of quarks carries a baryonic current which is same as 
the topological current. Thepossible connections of QCD to effective meson 
theories are due to 't Hooft^ and Witten^. They showed that in the limit 
of large number of colors iV c , QCD can be reduced to a theory of weakly 
interacting mesons and glue balls. Using 1/N C as an expansion parameter, 
Witten argued that to leading order, baryons emerge as solitons in the weakly 
coupled phase of mesons if some topological properties and stability condi- 
tions are assumed. Later Witten^ showed that it is necessary to introduce 
a Wess-Zumino term in the effective chiral theory in order to predict some 
observed physical processes allowed by the non-abelian anomaly in QCD. He 
also demonstrated that the gauged W-Z action can lead to the identifica- 
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tion of the topological current with the baryonic current. To verify another 
Skyrme conjecture about topological solitons being fermions, he rotated a 
skyrmion adiabatically through 27r-radians and found that the skyrmion had 
picked up a phase of N c tc. Hence a skyrmion is a fermion if N c is odd. 

Afterwards, Adkins, Nappi and WitteiJ^ systematically studied some 
static properties of nucleons in an explicit model in which the Skyrme La- 
grangian is invariant under SU(2) x SU(2). These results agreed with ex- 
periment to 30 %. Since then, there has been a tremendous interest in 
the Skyrme model. A good number of modifications to the model have 
been proposed, such as including the effect of the pion mass^^, adding 
higher order terms^l extending the symmetry to 5* £7(3) x SU(3) ^, in- 
cluding vector mesonsll^l, hidden symmetry approach bound state ap- 
proach to strangeness^^ , and non linear supersymmetric models and relation 
11 7111 ^1 

to skyrmions 1 11 *. A recent review of supersymmetric chiral models was 
given by PerelomovEH. 

Although the Skyrme model provides a reasonably good description of 
baryons, there are many problems which remain to be solved. In the Skyrme 
model, if the masses of the low lying baryons N and A are taken as input 
parameters, the predicted value of the pion decay constant F n is found to be 
too high; or if the experimental value of F w is taken as an input, the predicted 
masses of the nucleons turn out to be too high. Because the value F^ is too 
large, many other predictions which depend on F v disagree with experiment. 
The model also predicts a tower of baryon states with equal spin and isospin 
(s = i = 5/2,7/2,...) which do not exist in nature. They are considered 
artifacts of the model. 

Quark model, on the other hand, is regarded as an effective theory when 
the gluon degrees of freedom are integrated out. If both the Skyrme model 
and the quark model are effective models of QCD, they should have common 
symmetry of hadrons, such as chiral symmetr y T he hadronic approximate 
SU(6) symmetryEHl and SU (6/21) symmetryETJ have been established for 
the quark model but ways to restore these symmetries for the Skyrme model 
have yet to be found. 

In this paper we address some of these problems, give partial solutions 
to some of the difficulties by improving the Skyrmion model for hadrons 
and suggest a strategy toward building a full fledged quantitatively viable 
skyrmionic description of hadronic physics. 
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1 Effective Hadronic Super symmetry 



Supersymmetry is a symmetry between fermions of half odd integer spin and 
bosons of integer spin. In supergravity, for example, we have an action that 
is invariant under operations transforming fields of spin 0, |, 1, |, and 2 

among themselves^! . Here we recall that low lying hadrons have the same 
range of spin values with s = \ and § baryons interlaced with s = 0, 1 and 
2 mesons. Groups that transform s = and ~ mesons into each other and 

simultaneously mix s — \ and | baryons have been proposed long ago-"^-. 
This kind of symmetry had to be broken symmetry since no degeneracy 
is observed between particles of different spin. However, a cursory glance 
at the Chew-Frautschi plot of hadronic trajectories^^ reveals the following 
features: 

• (a) All trajectories with j (hadronic spin) versus m 2 (m = hadronic 
mass) are approximately linear. 

• (b) Leading mesonic trajectories associated with lowest spin 0, 1 and 2 
are parallel among themselves. 

• (c) Leading baryonic trajectories associated with s = ~, | are also 
parallel among themselves. 

• (d) Mesonic and baryonic trajectories are approximately parallel to 
each other with a universal slope of a ~ 0.9 (GeV)~ 2 . 

• (e) The separation between mesonic trajectories is nearly the same as 
the one between baryonic trajectories. 

Now, properties (b) and (c) suggest the existence of a phenomenological 
symmetry between mesons of different spin, which also operates on baryons 
with different spin. For hadrons composed of light quarks u, d and s, this 
symmetry is expressed by the group SU(6) x 0(3) where 0(3) describes the 
rotational excitations on the leading trajectories, and the spin-flavor group 
SU(6) classifies the lowest elements of the trajectories into particle multi- 
plets. The property (d) on the other hand, tells us that there is a new kind 
of symmetry (supersymmetry) between the bosonic mesons and fermionic 
baryons. The universal Regge slope for hadrons is a supersymmetric observ- 
able. 
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The meaning of property (e) is that, the physical mechanism that breaks 
the SU(6) symmetry must also be responsible for breaking of its supersym- 
metric extension. Finally, from property (a) we infer that the potential 
binding the quarks is approximately linear and we have to apply relativistic 
quantum mechanics appropriate to light quarks while the Schroedinger non- 
relativistic theory is sufficient for the description of quarkonium systems for 
heavy quarks. 

According to the quark model of Gell-Mann and Zweig, mesons and 
baryons are described respectively by bound (qq) and (qqq) systems. Any 
symmetry between mesons and baryons must correspond at the quark level to 
a supersymmetry between q (antiquarks) and bound (qq) states (diquarks). 
Now q, with s = | and unitary spin associated with the triplet (3) represen- 
tation of the flavor 577(3) belongs to the (6) representation of SU(6). The 
low lying baryons are in its (56) representation. Since (56) is contained in 
6 x 21 = 56 + 70, the diquark with s = or 1 must be in the symmetric (21) 
representation of SU(6). Hence the hadronic supersymmetry we are seeking 
must transform the (6) and (21) SU(6) multiplets, both color antitriplets 
into each other and therefore must be 27 dimensional with 6 fermionic and 
21 bosonic states. This supergroup is SU(6/21). It was first introduced by 
Miyazawa^^ as a generalization of the hadronic SU (6) symmetry, following 
earlier attempts by Hwa and Nuyts-^-. 

The (q) — (qq) symmetry that also implies (q) — (qq) symmetry will trans- 
form the meson (qq) in general not only to baryons (qqq) and antibaryons 
(qqq) but also to exotic mesons D — D (D = qq, and D = qq) that belong 
to the SU(6) representations 1,35 and 405. The 1 and (35) are + and 1 + 
mesons while the (405) also includes mesons with spin 2 + and isospin 2. All 
the low lying hadrons will now be in the adjoint representation of SU(6/21) 
with both spin and isospin taking values 0, |, 1, 1 and 2. 

The next introduction of supersymmetry into physics was within the con- 
text of dual resonance models that evolved later into string models. These 
theories due to Ramond^l and Neveu and Schwarz^l lead naturally to lin- 
ear baryon and meson trajectories that are parallel. The string models were 
not local and they were not relativistic in four dimensions. 

Examples of renormalizable relativistic local quantum field theories in- 
volving s = 0, \ and 1 fields with interaction were first constructed by Wess 

and Zumino ^^- following initial attempts by other authors ^ . Wess and 
Zumino based their work on the super-Poincare algebra that is a supersym- 
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metric generalization of the infinitesimal Poincare group and sits in the super- 
conformal algebra. The final step of incorporating s = 2 and s = | fields in a 

local relativistic field theory was taken by the discoverers of supergravity^l 
and extended supergravity who were able to supersymmetrize Einstein's gen- 
eral relativity and Kaluza-Klein theories. All these local supersymmetric 
quantum field theories have fascinating convergence and symmetry proper- 
ties but they are far from describing properties of hadrons or even quarks 
and leptons. If these fundamental fields have any physical reality at all, 
they may be associated with preons (or haplons) that would be hypothetical 
constituents of quarks, leptons and fundamental gauge bosons. 

Returning to the more concrete world of hadrons, we may try to see if the 
phenomenological approximate group SU(6) and its supersymmetric exten- 
sion can be justified within the standard theory of colored quarks interacting 
through gluons that are associated with the color gauge group SU(3) C . The 
justification of SU (6) and the derivation of its breaking was given by Georgi, 
Glashow and de Rujula 1 ^. 

On the other hand, a phenomenological supersymmetry in nuclear physics 
between odd and even nuclei was discovered and formulated by means of 
supergroups'^ U{m/n) by Balantekin, Bars and Iachello'"^ . It was natural 
to see if such groups could also describe hadronic supersymmetry. In the first 
part of this paper we propose to extend the approach of Georgi et al. to see 
how far it can provide a basis for the existence of an approximate hadronic 
supersymmetry. It turns out that most of the ingredients for this approach 
are already in the literature. 

The key concepts are: 

• (a) The approximate validity of the string theory as an approx- 
imation to QCD. This was shown by the lattice gauge theory as a 
strong coupling approximation to QCD, following the pioneering work 
of Wilson^S] anc l a i so by the elongated bag model of Johnson and 



ThornESl following the bag model" 3 ^ approximation to QCD. 

(b) The emergence of the diquark structure in QCD through the 
string approximation. This was done by Eguchi^^ and also by Johnson 
and Thorn 

(c) The vector (spin dependent) nature of the Coulomb part together 
with the scalar (spin independent) nature of the confining part of the 
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q — q potential. These properties were collectively worked out by 
many authors using both perturbation theory and lattice gauge the- 
ory methodsESl. 

• (d) The necessity for introducing exotic D — D mesons in QCD through 
the bag approximation^^, the string picture^^'E^ or the confining 
potential model^T 

• (e) Deviations from linearity of Regge trajectories within the context 

of qcdE3I. 

Once all the parts of this jigsaw puzzle is put together, a rather sim- 
ple picture emerges. The diquark behaves very much like the antiquark in 
the strong coupling regime because QCD forces are flavor independent and 
the confining part of QCD potential is spin independent. This immediately 
leads to a q — (qq) effective supersymmetry at large separation. At short 
distances, spin is approximately conserved because of asymptotic freedom 
but spin independence gets broken through one-gluon exchange in the per- 
turbation theory regime. The diquark structure also disappears at short 
distance, leading to the breaking of both 577(6) and SU (6/21) at low ener- 
gies. There is another difficulty associated with a supersymmetric extension 
of SU(6) noted by Salam and Sthrathdee^l. Because of the anticommuta- 
tivity of the Grassmann numbers and also of supercharges, antisymmetrical 
representations of SU(6) like (15) for diquarks and (20) for baryons will 
occur in hadronic supermultiplets. In the case of colored quarks, however, 
it is possible to introduce non associative Grassmann numbers Ui,U2 and 
constructed out of octonion units^^. They transform like a triplet under 
the SU(3) subgroup of the automorphism group G<i of octonions. Then the 
octonionic quarks u ■ q a and u ■ will commute unlike q l a and q^ that are 
anticommutative (a and (3 are combined spin-flavor indices). This procedure 
which is the basis of color algebra not only converts antisymmetrical spin- 
flavor group representations into symmetrical ones, but suppresses the color 
sextet (6) representations for diquarks that would otherwise be allowed un- 
der the rule for the preservation of the Pauli principle for the colored quark 
states. With the introduction of octonionic quarks, diquarks, mesons and 
baryons can all be viewed as elements of an octonionic superalgebraEB. 

The simplest supersymmetric Hamiltonian is obtained starting from semirel- 
ativistic dynamical models of quarks and diquarks already used by Lichten- 
berg et al.,^J for an approximate calculation of baryonic masses. In the 
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following subsections we will first give a derivation of hadronic mass for- 
mulae using semi-relativistic and relativistic formulation. Later, in section 
2 using the skyrmion model we will show that we obtain the same mass 
formulae. 

The Diquark and the Semirelativistic Hamil- 
tonians for (qq) and q — (qq) Systems 

In QCD, both q — q and q — q forces are attractive, the color factor of the 
latter being half of that of the former. Thus, the formation of diquarks inside 
a baryon is a definite possibility. Various strong coupling approximations to 
QCD, like lattice gauge theoryE^ \ Hooft's approximation (when N c , 

the number of colors is very large), or the elongated bag model^^ all give 
a linear potential between widely separated quarks and an effective string 
that approximates the gluon flux tube. In such a theory, EguchiEU has 
shown that it is energetically favorable for three quarks in a baryon to form 
a linear structure with a quark in the middle and two at the ends or, for high 
rotational excitation, a bilocal linear structure with one quark at one end 
and a diquark at the other end. This was reconfirmed by Johnson and Thorn 
independently in the bag model when the bag is deformed and elongated in 
a rotationally excited baryon Thus if we move along a leading baryon 
Regge trajectory to a region of high j, we are likely to find a baryon as 
a bilocal object consisting of a quark and a diquark instead of a trilocal 
object that represents a ground state baryon more accurately. On the other 
hand a meson is a bilocal object consisting of a quark and an antiquark 
interacting via a linear potential at large separation. Consider then a bilocal 
object with constituents with respective masses mi, m 2 , spins s 1; s 2 and color 
representations 3 and 3. One of the constituents can be a quark q or the 
anti-diquark D (both color triplets) while the other can be an antiquark q or 
a diquark D which are both color antitriplets. The QCD force between the 
two constituents will be flavor independent and it will consist of two parts: 
a Coulomb like part V c which transforms like the time component of a 4- 
vector (due to the exchange of single s = 1 gluons at short separations) and 
a confining part V s which is largely a relativistic scalar, is spin independent 
and is due to the exchange of a great many gluons that form a flux tube or 
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an elongated bag at large separation of the constituents. 

Let Pi = p and p 2 = —p be the center of mass momenta of the two 
constituents. The quantity which is canonically conjugate to the relative 
coordinate r = r± — r 2 is 

d 1 

~ l v = ~ l = 2 ~ P2) =p ^ 

Ignoring the center of mass motion, following Lichtenberg et al.^^, we can 
write a semi-relativistic wave equation for the wave function i\)\ 2 (r) of the 
bilocal object with energy eigenvalues Qi 2 , namely 



(Q 12 - V c )4j 12 = {[(m a + ± V s ) 2 - V 2 ]U [(m 2 + ± K) 2 - V 2 )]^ 12 (2) 
The scalar and vector potentials are given by 

V s = br, V c = h Kia (3) 

where | is the color factor, ot s is the strong coupling constant at the energy 
fii2, and the spin dependent part of the vector potential is the hyperfine 
structure correction due to gluon exchange with k 12 = 1^12 (0) 1 2 - We see 
that at large r, neglecting the mass difference (m 2 — mi), we find the same 
equation for both the (q — q) and the q — D system, except for the presence 
of the hyperfine term that breaks the symmetry between q and D. To this 
approximation, we can transform the second constituent q into D and vice- 
versa without changing the energy eigenvalue Q. This means that the system 
admits the approximate SU (6/21) supersymmetry transformation 

S& = b aPl (D^ 

5(D i ) l3y = b a ^ qi (4) 
in addition to the SU(6) transformation 

8{D i )^ = n^l (Dy* (5) 
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The breaking of both SU(6) and SU (6/21) is due to the hyperfine term 
while the supersymmetry is further broken by the quark-diquark mass dif- 
ference mi — TO 2 . 

We could also have brought into play the wave functions ipo and tp q of 
the diquark and antidiquark at point r 2 in the field of the quark at point r 2 . 
The masses mi and m% must then be replaced by the reduced masses 

m q m D m q m q 1 

= : , fiq = : = - m q (6) 

m q + m& rn q + m q 2 

In this case, the wave function belongs to the fundamental 27 dimen- 
sional representation of SU(6/21) and the hamiltonian commutes with the 
supersymmetry transformation Eq.Q of the wave function except for the 
difference in the reduced masses \i£> and \i q and the hyperfine structure term 
in the Hamiltonian. 

The breaking of supersymmetry results in the mass difference 

Am = m 2 — mi + k — (7) 

mim2 

at high energies both for baryons and mesons. At low energies, the baryon 
mass becomes a trilocal object with three quarks and the mass splitting is 
given by 

1 , Si ' So So ' Si • Si . , . 

Ami23 = 77 K + + 8 

2 mim 2 m 2 m3 m^mi 

where mi, m 2 and m^ are the masses of the three different quark constituents. 

Going back to the first formulation which is more symmetrical, we can 
write the Hamiltonian associated with Eq.(j2J) in a SU(6/21) covariant form: 
we have four cases for t/> 12 , namely ip qq , ip^n, ip q u and ipQ q that represents 
mesons, exotic mesons, baryons and antibaryons with respective dimensions 
6 x 6, 21 x 21, 6 x 21, and 21 x 6 that all fit into the following 27 x 27 adjoint 
representation of SU(6/21): 

^ - ( t& ) ^ 

The wave equation for the hadronic wavefunction ip can now be written 

as 

i^ = Hil> = [K,4) + Sil>S (10) 
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where K is the diagonal matrix 



K = ( ^ q ^ ^j(2i) ) 7 1^ = nx n unitmatrix (11) 



with 



K D = - 2 -^ + [(m q+ l -V s f-W^ (12) 



and 



The second, spin dependent term on the right hand side of Eq. ()l()jl is 
a symmetry breaking term for both SU(6) and SU(6/21). The first term 
is SU(6) symmetrical. It also preserves the supersymmetry in the limit 
m q = m D . Hence if the quark-di quark mass difference and the spin dependent 
terms are neglected, Eq.(|l(J|) is invariant under the SU(6/21) infinitesimal 
transformation 

8^=[N,ip] (14) 

where 

"-(I'm) < 15 > 

where M and N have elements mP a and n^J respectively in Eq.jjHJ) while 
the rectangular matrices b and b have elements b a ^ and 6 Q/ g 7 occurring in 
Eq.(j3|. In that limit, the Hamiltonian H commutes with the generators of 
the transformation Eq. (jl4j) . 



The Hadronic Regge Trajectories 

Consider the Hamiltonian of Eq. (J2J) . We can write 



V =p =p r + — 2 (16) 
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where i is associated with the orbital excitation of the system. For high 
rotational excitations, the expectation value of r is large, corresponding to 
a stretched string. The angular momentum i is also large. The value of the 
centrifugal energy which is proportional to has a similarly large value. 

Since V s that is proportional to r will also have a high absolute value, the 
constituent masses become negligible in the high relativistic limit. On the 
other hand, the radial excitation term p 2 . can be neglected on the leading 
trajectory associated with the lowest radial energy. 

The ground state energy eigenvalue E of the Hamiltonian can be es- 
timated by using the Heisenberg uncertainty principle. This leads to the 
replacement of r by Ar and p r by 

A Pr = i(Ar)-\ (h=l). (17) 

Then E as a function of Ar is minimized for the value of r of Ar. 
The r corresponds to the Bohr radius for the bound state. The confining 
energy associated with this Bohr radius is obtained from the linear confining 
potential S(r) = br, so that the effective masses of the constituents become 

Mi =7711 + ^0, M 2 = m 2 + ^S , {S = br ) (18) 

For a meson mi and m 2 are the current quark masses while Mi and M 2 can 
be interpreted as the constituent quark masses. Note that even in the case of 
vanishing quark masses associated with perfect chiral symmetry, confinement 
results in non zero constituent masses that spontaneously break the SU (2) x 
SU (2) symmetry of the u, d quarks. 

Let us illustrate this method on the simplified spin free Hamiltonian in- 
volving only the scalar potential. In the center of mass system, = 0, 
or = — p^ = p. The semi-relativistic hamiltonian of the system is then 
given by 

E 12 <S> = J2^(m l+ l -brf + p 2 $. (19) 
Taking m x = m 2 = m for the quark- ant iquark system, we have 

E l2 $ = 2^(m+i&r)2 + p2 + ^(^±i) ( 20 ) 
where we have written the momentum part in spherical coordinates. 
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Putting 

b = p 2 , p = p r, (21) 
for the q — q system we find Eyi by minimizing the function 



E qq =2 x l{m+ l -^f + ^+ l -Y. (22) 



For u and d quarks, m is small and can be neglected so that 

E 2 = p 2 [p 2 + p- 2 (2£ + l) 2 } (23) 
which has a minimum for 

p 2 = p 2 = 2£+l, (24) 

giving 

E 2 mm = E 2 {p Q )=Ap 2 {i+ 1 -). (25) 

Thus, we obtain a linear Regge trajectory with 

'lob 
a = -p~ 2 = -. (26) 

Also J = I + S , where S arises from the quark spins. Experimentally 

a = 0.88(GeV)~ 2 (27) 

for mesons giving the value 0.54 GeV for p. A more accurate calculation (see 
[4]) gives 

a = (2vr/i 2 )- 1 , p ~ OAZGeV. (28) 
The constituent quark mass can be defined in two ways 



M c (e) = ^E min = pJe + ^, (29) 



or 



m' c (£) = S = ^p Po = ^ i + X -. (30) 
The first definition gives for £ = 0, 

M c = 0.3lGeV for ^ = 0.43 (31) 
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in the case of u and d quarks. 

When the Coulomb like terms are introduced in the simplified Hamilto- 
nian 1)20)1 with negligible quark masses one obtains 

E = £ [-a + v / P 4 + (2£+l) 2 ] (32) 

with 

4 2 
a = -a s for (qq), a = -a s for (qq). (33) 

In the energy range around 1 GeV, a s is of order of unity. Estimates range 
from 0.3 to 3. Minimization of E gives 

E = fxuf (-a + v / % + (2£+l) 2 ) (34) 

where 

«„(£) = pt= (21 + 1) 2 (1 + \f + eV2/3)Je + yP), (35) 

e = ±1 - i3 =(2?TT)- (36) 

The minimum Eq is obtained for e = — 1, giving to second order in f3: 
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E = ^2(21 + 1)(1 - ^= - 3^). (37) 

Linear Regge trajectories are obtained if f3 2 is negligible. Then for mesons 

El = A/j, 2 £ + 2/i 2 (l - V2a). (38) 

The f3 2 is negligible for small £ only if we take the lowest estimate for a s , giv- 
ing 0.4 for a in the qq case. For mesons with u, d constituents, incorporating 
their spins through the Breit term we obtain approximately 

c 3c AV 

m p ~m w = £ +-, m n = E - — , c = K — (39) 

where M g is the constituent quark mass. This gives 

E = (3m ^ +m - } = OMGeV. (40) 
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The Regge slope being of the order of lGeV an average meson mass of the 
same order is obtained from Eq. (J37|) in the linear trajectory approximation. 
To this approximation a should be treated like a parameter rather than be 
placed by its value derived from QCD under varying assumptions. Using 
Eq. (j2~%|) for /i one gets a better fit to the meson masses by taking a s ~ 0.2. 

Turning now to baryon masses, we must first estimate the diquark mass. 
We have for the qq system 

M D = /i( v / 2-^« s ), (41) 
that is slightly higher than the average meson mass 

m = fx(V2-^a s ). (42) 

Here we note that E is not very sensitive to the precise value of the QCD 
running coupling constant in the GeV range. Taking a s ~ 0.3 changes E qq 
from 0.55 to 0.56GW. 

Note that Eq. (j4*T|) gives m D = 0.55GeV. For excited q — q and q — D 
systems if the rotational excitation energy is large compared with /i, then 
both the rap and the Coulomb term — (same for q — D and q — q systems) 
can be neglected. Thus, for both (q — D) [excited baryon] and q — q [excited 
meson] systems we have Eq . (J25j) . namely 

(E q - D f ~ (E q ~ q f ~ 4/i 2 £ + 2/i 2 (43) 

giving again Eq. (}2T)j) . i.e. 

1 

47T 2 " V 27T/i 2 



(«V»=(«')h=^ or (^) ( 44 ) 



as an explanation of hadronic supersymmetry in the nucleon and meson 
Regge spectra. We also have, extrapolating to small £: 

A(M 2 ) q ~ D = A(m 2 ) q - q = 4fi 2 Ai = ^A£. (45) 

a 



m A ~ m N = m p — m l- (46) 



For Ai = 1 we find 



This relationship is same as the one proved in our earlier paper 1 ^ through 
the assumption that U (6/21) symmetry is broken by an operator that behaves 
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like s = 0, / = member of 35 x 35 representations of SU(6), which is true 
to 5%. It corresponds to a confined quark approximation with a s = 0. 

The potential model gives a more accurate symmetry breaking (a s ~ 0.2): 

9 

-{m 2 p - ml) =m\-m 2 N (47) 

with an accuracy of 1%. 

This mass squared formula arises from the second order iteration of the 
q — D, q — q Dirac equation. The factor | comes from 

Vl as)2 = \* (48) 

At this point it is more instructive to derive a first order mass formula. 
Since the constituent quark mass M q is given by Eq.(|2HJ) (£ = 0), we have 

M, = £ (49) 

so that 

\/2 

m = 2M q (l - —a s ) ~ 1.9M q (50) 
3 

When the baryon is regarded as a q — D system, each constituent gains an 
effective mass \[ipo which was approximately the effective mass of the quark 
in the meson. Hence, the effective masses of q and D in the baryon are 

m q ~ M q , m D = M D + M q ~ 3M q (51) 

The spin splittings for the nucleon N and the A are given by the Breit 
term 

AM = KAV^^ (52) 
m q m D 

For the nucleon with spin | the term S^-S^ gives —1 while it has the value 
| for A with spin |. Using the same K for mesons and baryons which are 
both considered to be a bound state of a color triplet with a color antitriplet 
we can relate the baryon splitting AM to the meson splitting Am for which 
S q ■ S q takes the values \ and Hence we find 

aw w w 3 KAV 1 ^ AV " i a KAV /roN 

AM = M A -M N = -.^ = -.^ | -, and Am = (53) 



16 



which leads to a linear mass formula 



AM = ^Am (54) 

which is well satisfied, and has been verified before using the three quark 
constituents for the baryonE^. 

The formation of diquarks which behave like antiquarks as far as QCD is 
concerned is crucial to hadronic supersymmetry and to quark dynamics for 
excited hadrons. The splittings in the mass spectrum are well understood 
on the basis of spin-dependent terms derived from QCD. This approach to 
hadronic physics has led to many in depth investigations recently. For exten- 
sive references we refer to recent papers by Lichtenberg and collaborators^^ 
and by Klempt^^. To see the symmetry breaking effect, note that the mass 
of a hadron will take the approximate form 

m 12 =mi + m 2 + K— (55) 

mi77l2 

where rrii and Sj (i = 1, 2) are respectively the constituent mass and the spin 
of a quark or a diquark. The spin-dependent Breit term will split the masses 
of hadrons of different spin values. If we assume m q = m 9 - = m, where m is 
the constituent mass of u or d quarks, and denote the mass of a diquark as 
mo, then this approximation gives 

3 1 

= (m q q) s = = 2m - K ^2i m p = (w%) s =i = 2m + K^—^ (56) 

m A = (m qD ) s=3/2 = m + m D + K- — — (57) 

ImrriD 

m N = {m qD ) s= i/ 2 = m + m D - K — - — , (58) 

mrrif) 

Eliminating m, nip and K, we obtain a mass relation 
8 2m A + m N 3 m p - m n 

o ' ^ , = 1 + o ^ 59 J 

3 3m p + m w 2 m A — m^ 

which agrees with experiment to 13%. 
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Relativist ic Case 

Let us consider a quark-antiquark system in the approximation that the 
potential is only a scalar. In the center of mass system, p^ + p^ = 0, or 
p^ = —p^ = p The semi-relativistic hamiltonian of the system is then 
given by 



i=l 



ff = EvK + l br ) 2 + p 2 ( 60 ) 



Taking m x = m 2 = m for the quark-antiquark system, we have 



H= 2\f {m + -br) 2 + p 2 r + v - ' (61) 



where we have written the momentum part in spherical coordinates. For u 
and d quarks, m is small and can be neglected so that 



F 2 = 4[ l 6 2 r 2 +p 2 + ^+l) ] (62) 



Similar equations for aq — q meson system were already proposed and solved 
numerically^^ ' or exactly^! by several authors. The eigenfunction \I/ for 
H 2 has eigenvalue E 2 such that 

r i, 2 2 i d 2 e(£ + i), T ^ 2t 

The differential equation can be solved exactly, and and the normalized eigen- 
function is found to be 



2(f) 2 ^ 3 r(£ + § + n r -l) 



l)\T 2 {i + 



1/2 



f r u 2n 

r exp| — — r J x 



F(-n r + l,£ + Hr 2 )Yr( 



(64) 



where F(—n r + 1,£+ |, |r ) is the confluent hypergeometric function, and 
n r = 1, 2, 3, . . ., is the radial quantum number. The eigenvalue is given by 



E 2 = Ab[2(n r 



(65) 
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therefore we obtain linear Regge trajectories of slope 1/46 when we make 
plots of £ versus M 2 . The case of n r = 1 corresponds to the leading Regge 
trajectory, and cases of n r = 2,3,..., correspond to the parallel daughter 
trajectories. 

Note that Iachello and his collaborators^^ have obtained a similar mass 
formula based on algebraic methods. Starting from a spectrum generating 
algebra G, we can write a chain of subalgebras 

GDG' DG" D ■■■ . (66) 

They proposed that the Hamiltonian can be expanded in terms of invariants 
of the chain of subalgebras 

H = aC{G) + f3C(G') + 7 C(G") + ■ ■ • , (67) 

where C(G) denotes one of the invariants of G. In a two-body problem, the 
spectrum generating algebra is Z/(4), and in the total symmetric represen- 
tation, one of the two chains of subalgebras is U(3) followed by SO (3) with 
respective Casimir invariants n(n + 2) and l{l + 1), where n corresponds to 
the vibrational mode. In order to derive the relativistic mass formula like 
Eq.(|65|). one has to write a formula for the square of the mass, rather than 
the mass itself. Furthermore, without violating the dynamic symmetry, the 
Hamiltonian can be written in terms of non-linear functions of the Casimir 
invariants. Therefore, we can have 



M 2 = a + Py/n(n + 2) + 1 + 7v /£(£ + 1) + 1/4 

= a' + pn + ji (68) 

Comparing this mass formula with Eq. (J65|) . we find that the relativistic quark 
model suggests a relation of (3 = 2 7 in Eq.(|55j). However, the above algebraic 
method treats (3 and 7 as two independent parameters, and so the extra 
degree of freedom allows Eq.()68|) a better fit to the experimental data. The 
fitted parameters for mesons are found to be: (3 ~ 1.5 (GeV) 2 and 7 ~ 1.1 
(GeV) 2 . 

The inclusion of the Coulomb-like term causes deviations from the linear 
trajectories at low energies and changes the relation between the parameters 
(3 and 7 . It is important to carry out a calculation based on a better ap- 
proximation to see if the empirical values of these parameters are compatible 
with our model. 
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Bilocal Approximation to Hadronic Structure 
and Inclusion of Color 



The supergroup SU (6/21) acts on a quark and antidiquark situated at the 
same point T\. At the point r 2 we can consider the action of this supergroup 
with the same parameters, or one with different parameters. In the first case 
we have a global symmetry. In the second case, if we only deal with bilocal 
fields the symmetry will be represented by SU (6/21) x SU(6/21), doubling 
the supergroup. On the other hand, if any number of points are considered, 
with different parameters attached to each point, we are led to introduce a 
local supersymmetry SU (6/21) to which we should add the local color group 
SU(3) C . Since it is not a fundamental symmetry, we shall not deal with the 
local SU (6/21) group here. However, the doubling the supergroup is useful 
for bilocal fields since the decomposition of the adjoint representation of the 
728-dimensional SU "(6/21) group with respect to SU(6) x SU (21) gives 

728 = (35, 1) + (1, 440) + (6, 21) + (6, 21) + (1, 1) (69) 

A further decomposition of the double supergroup into its field with re- 
spect to its center of mass coordinates, as will be seen below leads to the de- 
composition of the 126-dimensional cosets (6, 21) and (21, 6) into 56 + + 70~ 
of the diagonal SU(6). We would have a much tighter and more elegant 
scheme if we could perform such a decomposition from the start and be able 
to identify (1,21) part of the fundametal representation of SU(6/21) with 
the 21-dimensional representation of the SU(6) subgroup, which means going 
beyond the SU (6/21) supersymmetry to a smaller supergroup having SU(6) 
as a subgroup. Full description of a bilocal treatment and a minimal scheme 
where the bilocal treatment gets carried over unchanged into it will be a 
subject of another publication. 

Setting up of an effective relativistic theory based on the Dirac equation 
for the quark and the Klein-Gordon equation for the antiquark, exhibiting 
invariance under relativistic supersymmetry is needed. This can be achieved 
by means of an effective Wess-Zumino type Lagrangian constructed out of 
a gluon field interacting with a vector superfield formed by a quark and an 
antidiquark. In what follows we will present some preliminary work in this 
direction. 

Low-lying baryons occur in the symmetric 56 representation^^ of SU(6), 
whereas the Pauli principle would have led to the antisymmetrical 20 rep- 



20 



resentation. This was a crucial fact for the introduction of color degree of 
freedom^! based on SU(3) C . Since the quark field transforms like a color 
triplet and the diquark like a color antitriplet under SU(3) C , the color degrees 
of freedom of the constituents must be included correctly in order to obtain 
a correct representation of the q-D system. Hadronic states must be color 
singlets. These are represented by bilocal operators 0(ri,r 2 ) in the bilocal 
approximation^S! that gives g(l)g(2) for mesons and -D(l)g(2) for baryons. 
Here g(l) represents the antiquark situated at r 1; q(2) the quark situated at 
r 2 , and D{\) = q{l)q{l) the diquark situated at ri. If we denote the cm. 
and the relative coordinates of the consituents by R and r, where r = r 2 — ri 
and 

R = , \ ' (70) 

(mi + m 2 ) 

with irti and m 2 being their masses, we can then write 0(R, r) for the op- 
erator that creates hadrons out of the vacuum. The matrix element of this 
operator between the vacuum and the hadronic state h will be of the form 

</i|O(R,r)|0>= X (RMr) (71) 

where x(R) is the free wave function of the hadron as a function of the cm. 
coordinate and ip(r) is the bound-state solution of the £7(6/21) invariant 
Hamiltonian describing the q — q mesons, q — D baryons, q — D antibaryons 
and D — D exotic mesons, given by 



iMap = [\ {m a + \v s y + p 2 + J (m p + \v s ) 2 + p 2 - ^ + fc^L^]^ 
V 2 V 2 3 r m a m/3 

(72) 

Here p = —iV in the cm. system and m and s denote the masses and 
spins of the constituents, a s the strong-coupling constant, V s = br is the 
scalar potential with r being the distance between the constituents in the 
bilocal object, and k = |^(0)| 2 . 

The operator product expansion^! will give a singular part depending 
only on r and proportional to the propagator of the field binding the two 
constituents. There will be a finite number of singular coefficients c n (r) 
depending on the dimensionality of the constituent fields. For example, for 
a meson, the singular term is proportional to the progagator of the gluon 
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field binding the two constituents. Once we subtract the singular part, the 
remaining part 0(R, r) is analytic in r and thus we can write 

0(R,r) = O (R) + r-Oi(R) + O(r 2 ). (73) 

Now Oo(R) creates a hadron at its cm. point R equivalent to a £ = 0, 
s-state of the two constituents. For a baryon this is a state associated with q 
and D ~ qq at the same point R, hence it is essentially a 3-quark state when 
the three quarks are at a common location. The Oi(R) can create three 
£ = 1 states with opposite parity to the state created by Oo(R). Hence, if 
we denote the nonsingular parts of g(l)g(2) and _D(l)g(2) by [g(l)g(2)] and 
[-D(l)g(2)], respectively, we have 

[g(l)g(2)]|0 >= |M(R) > +r • |M'(R) > +0(r 2 ) (74) 

[L>(l)g(2)]|0 >= \B(R) > +r • |B'(R) > +0(r 2 ) (75) 

and similarly for the exotic meson states D(1)D(2). 

Here M belongs to the (35 + 1) -dimensional representation of SU(6) cor- 
responding to an £ = bound state of the quark and the antiquark. The 
M (R) is an orbital excitation [£ = 1) of opposite parity, which are in the 
(35 + 1,3) representation of the group SU(6) x 0(3), 0(3) being associated 
with the relative angular momentum of the constituents. The M states con- 
tain mesons like B, A x , A 2 and scalar particles. On the whole, the £ = 
and £ = 1 part q(l)q(2) contain 4 x (35 + 1) = 144 meson states. 

Switching to the baryon states, the requirement of antisymmetry in color 
and symmetry in spin-flavor indices gives the (56) + representation for £>(R). 
The £ = 1 multiplets have negative parity and have mixed spin-flavor sym- 
metry. They belong to the representation (70~,3) of SU(6) x 0(3) and are 
represented by the states |B (R) > which are 210 in number. On the whole, 
these 266 states account for all the observed low-lying baryon states obtained 
from 56 + 3 x 70 = 266. A similar analysis can be carried out for the exotic 
meson states _D(1)D(2), where the diquark and the antidiquark can be bound 
in an £ = or £ = 1 state with opposite parities. 
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Color Algebra and Octonions 

The exact, unbroken color group SU(3) C is the backbone of the strong inter- 
action. It is worthwhile to understand its role in the diquark picture more 
clearly. 

Two of the colored quarks in the baryon combine into an anti-triplet 
3 x 3 = 3 + (6), 3 x 3 = 1 + (8). The (6) partner of the diquark and the 
(8) partner of the nucleon do not exist. In hadron dynamics the only color 
combinations to consider are 3 x 3 — > 3 and 3x3^1. These relations 
imply the existence of split octonion units through a representation of the 
Grassmann algebra {ui,Uj} = 0, i = 1,2,3. What is a bit strange is that 
operators Uj, unlike ordinary fermionic operators, are not associative. We also 
have h[ui, Uj] = u* k . The Jacobi identity does not hold since [ui, [uj, Uk]] = 
—ie-j 7^ 0, where e?, anticommute with Ui and u*. 

The behavior of various states under the color group are best seen if 
we use split octonion units defined by^ 

1 1 

u = -(l + ie 7 ), u* = -(l-ie 7 ) (76) 

uj = he-j + ie j+3 ), u* = -(e, - ie j+3 ), j = 1, 2, 3 (77) 

The automorphism group of the octonion algebra is the 14-parameter 
exceptional group G<i- The imaginary octonion units e a {a = 1, 7) fall into 
its 7-dimensional representation. 

Under the SU(3) C subgroup of G2 that leaves e-j invariant, uo and are 
singlets, while Uj and u* correspond, respectively, to the representations 3 
and 3. The multiplication table can now be written in a manifestly SU(3) C 
invariant manner (together with the complex conjugate equations): 



U 2 Q = U , U U* q = (78) 

uoUj = UjUQ = Uj, u^Uj = UjUo = (79) 

UiUj = -UjUi = e ijk u* k (80) 

UiUj = -5ijU (81) 



23 



where e^k is completely antisymmetric with = 1 for ijk = 123, 246, 435, 
651, 572, 714, 367. Here, one sees the virtue of octonion multiplication. If 
we consider the direct products 

C : 3 <g> 3 = 1 + 8 (82) 

G : 3 ® 3 = 3 + 6 (83) 

for SU (3) c , then these equations show that octonion multiplication gets rid 
of 8 in 3 ® 3, while it gets rid of 6 in 3 <g> 3. Combining Eq.® and Eq.flEU) 
we find 

(UiUj)u h = -e ijk u* (84) 

Thus the octonion product leaves only the color part in 3<S>3 and 3<8>3<8>3, 
so that it is a natural algebra for colored quarks. 

For convenience we now produce the following multiplication table for the 
split octonion units: 





u 


u* 


u k 


K 




u 





Uk 





u* 





ul 





K 


Uj 





Uj 




-5 jk u 


* 

u 1 


* 

.7 





-5 jk u* 





It is worth noting that U{ and u* behave like fermionic annihilation and 
creation operators: 

{u h uj} = {u*, u*} = 0, {m, u* k } = -Sij (85) 

For more recent reviews on octonions and nonassociative algebras we refer 
to papers by Okubo 22, and Baez^l. 

The quarks, being in the triplet representation of the color group SU (3) c , 
are represented by the local fields q l a (x), where % = 1,2,3 is the color in- 
dex and a the combined spin-flavor index. Antiquarks at point y are color 
antitriplets q l p(y). Consider the two-body systems 

C% = i&MqlM (86) 
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G% = gi(*i)gj(* 2 ) (87) 

so that C is either a color singlet or color octet, while G is a color antitriplet 
or a color sextet. Now C contains meson states that are color singlets and 
hence observable. The octet q — q state is confined and not observed as 
a scattering state. In the case of two-body G states, the antitriplets are 
diquarks which, inside a hadron can be combined with another triplet quark 
to give observable, color singlet, three-quark baryon states. The color sextet 
part of G can only combine with a third quark to give unobservable color 
octet and color decuplet three-quark states. Hence the hadron dynamics 
is such that the 8 part of C and the 6 part of G are suppressed. This can 
best be achieved by the use of the above octonion algebra'^. The dynamical 
suppression of the octet and sextet states in Eq.tjHHJ) and Eq.(jH7j) is , therefore, 
automatically achieved. The split octonion units can be contracted with color 
indices of triplet or antitriplet fields. For quarks and antiquarks we can define 
the "transverse" octonions (calling u and Uq longitidunal units) 

q a = Uiq l a = u • q Q , q p = ujgj = -u* • q> (88) 

We find 

?«(1)^(2) = u q a (l) ■ q/3 (2) (89) 

g Q (l)g^(2) = n*q Q (l)- q/3 (2) (90) 

G aP {12) = q a {l)qp{2) = u* ■ q a (l) x c^{2) (91) 

G? /9a (21) = qp{2)q a {l) = u* ■ Q/3 (2) x q a (l) (92) 
Because of the anticomutativity of the quark fields, we have 

G aP (12) = Gp a {21) = i{g (l), q p {2)} (93) 

If the diquark forms a bound state represented by a field D a p{x) at the 
center-of-mass location x 

1 

x = -(x 1 + x 2 ) (94) 
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when %2 tends to x\ we can replace the argument by x, and we obtain 



D a p(x) = Dp a {x) 



(95) 



so that the local diquark field must be in a symmetric representation of the 
spin-flavor group. If the latter is taken to be 577(6), then D a p{x) is in the 
21-dimensional symmetric representation, given by 



If we denote the antisymmetric 15 representation by A a p, we see that 
the octonionic fields single out the 21 diquark representation at the expense 
of A a p. We note that without this color algebra supersymmetry would give 
antisymmetric configurations as noted by Salam and Strathdee^l in their 
possible supersymmetric generalization of hadronic supersymmetry. Using 
the nonsingular part of the operator product expansion we can write 



The fields A a/3 have opposite parity to D a/3 ; r is the relative coordinate at 
time t if we take t = t\ = t%. They play no role in the excited baryon 
which becomes a bilocal system with the 21- dimensional diquark as one of 
its constituents. 

Now consider a three-quark system at time t. The cm. and relative 
coordinates are 



(6® 6), = 21 



(96) 




(97) 




(98) 



1 



(2r 3 - n - r 2 ) 



(99) 



P 




(100) 



giving 




~FP + ~1= Y 



1 1 



(101) 



~ V3 V / 6 P 



(102) 
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The baryon state must be a color singlet, symmetric in the three pairs 
(a, xi), ((3, x 2 ), (7, x 3 ). We find 



(g Q (l) g/3 (2))g 7 (3) = -u* F af3j (123) 



(104) 



? 7 (3)(?«(l)«i9(2)) 



M F Q/37 (123) 



(105) 



so that 



-{{g a (l),g /3 (2)} ) g 7 (3)} = F Q/37 (123) 



(106) 



The operator F Q/37 (123) is a color singlet and is symmetrical in the three 
pairs of coordinates. We have 



where C is of order two and higher in p and r. Because R is symmetric 
in ri, r 2 and r 3 , the operator B a ^ that creates a baryon at R is totally 
symmetrical in its flavor-spin indices. In the SU (6) scheme it belongs to the 
(56) representation. In the bilocal q — D approximation we have r = so that 
F Q/ g 7 is a function only of R and p which are both symmetrical in and r 2 . As 
before, B' belongs to the orbitally excited 70~ represenation of SU(6). The 
totally antisymmetrical (20) is absent in the bilocal approximation. It would 
only appear in the trilocal treatment that would involve the 15-dimensional 
diquarks. Hence, if we use local fields, any product of two octonionic quark 
fields gives a (21) diquark 



and any nonassociative combination of three quarks, or a diquark and a quark 
at the same point give a baryon in the 56 + representation: 



F Q/37 (123) = B af3y (R) + p ■ B'(R) + r ■ B"(R) + C 



(107) 





(g Q (R) g/9 (R))g 7 (R) 



^0-Ba/3 7 (R) 



(109) 



g Q (R)( g/9 (R)g 7 (R)) 



(110) 
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g 7 (R)(g a (R) g/3 (R)) 



(111) 



(g 7 (R)g Q (R))g^(R) 



(112) 



The bilocal approximation gives the (35 + 1) mesons and the 70 ~ baryons 
with £ = 1 orbital excitation. 

In order to go beyond the SU (6/21) symmetry we could go to a smaller 
supergroup and build a colored supersymmetry scheme based on SU(6) x 
577(6/1), for example. Applications of such extended algebraic structures 
and their implication for the construction of supersymmetric meson-baryon 
lagrangians will be deferred to another publication. 

We can now consider the 28 x 28 octonionic matrix 



Here e can be given values 1, —1 or 0. M and N are respectively 6x6 
and 21 x 21 hermitian matrices, B a rectangular 6 x 21 matrix, u ■ Q a 
6x1 column matrix, u ■ D* a 21 x 1 column matrix, and L a 1 x 1 scalar. 
Such matrices close under anticommutator operations for e = 1. Matrices iZ 
close under commutator operations. In either case, they do not satisfy the 
Jacobi identity. But for e = 0, when the algebra is no longer semi-simple, 
the Jacobi identity is satisfied and we obtain hadronic superalgebra which 
is an extension of the algebra SU(6/21). Its automorphism group includes 
SU(6) x SU(21) x SU(3) C . Thus color is automatically incorporated. 



Relativistic Formulation Through the Spin Re- 
alization of the Wess-Zumino Algebra 



It is possible to use a spin representation of the Wess-Zumino algebra to 
write first order relativistic equations for quarks and diquarks that are in- 
variant under supersymmetry transformations. In this section we briefly deal 
with such Dirac-like supersymmetric equations and with a short discussion 
of experimental possibilities for the observation of the diquark structure and 




u M u B u ■ Q 
u B^ u N u ■ D* 
eu* ■ Q t eu* ■ D T u L 




(113) 
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exotic D — D = (qq){qq) mesons. For an extensive discussion of the exper- 
imental situation we refer to recent papers by Anselmino, et al.S3 ; and by 
KlemptHS. 

There is a spin realization of the Wess-Zumino super-Poincare algebra 
[p M ,p„] = 0, [D a ,p M ] = (114) 

[D $ ,pA = 0, = <t*V (115) 

with p^ transforming like a 4- vector and D a , D@ like the left and right handed 
spinors under the Lorentz group with generators J M „. 
We also note that 

[Jftv,P\] = &m>P\ ~ S v\Pn ( 116 ) 

and 

[J,J} = J (117) 

The finite non unitary spin realization is in terms of 4 x 4 matrices for 
.1^ and p v 

•V = i<V = ^[7m>7v] (H8) 

J%, = 1 ~Y^l^ = K ( 119 ) 



P , = nj = (120) 



Introducing two Grassmann numbers 9 a (a = 1, 2) that transforms like 
the components of a left handed spinor and commute with the Dirac matrices 
7 M , we have the representation 

D a = A a = J- (121) 
= A? = 9 a afu^ (122) 
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Such a representation of the super-Poincare algebra acts on a Majorana 
chiral superfield 



S(x, 6) = ^(x) + 6 a B a (x) + -ej a X (x). 



(123) 



Here ip and x are Majorana superfields associated with fermions and B a 
has an unwritten Majorana index and a chiral spinor index a, so that it 
represents a boson. 

Note that the sum of the two representations we wrote down is also a 
realization of the Wess-Zumino algebra. 

On the other hand we have the realization of in terms of the differential 
operator — id^ = — i J^. In the Majorana representation, the operator 7^9^ = 
ij^Pfj, is real, and if) = C = ijj* . Let us now define ipL and ipR by 



fa = ^ + 7s)0 



(124) 



and 



*pR = -(1-75)^ = Vl 



(125) 



The free particle Dirac equation can now be written as 




(126) 



or 



n M p M 0L = -imi/j* L 
We can introduce left and right handed component fields 



(127) 



S L = -(l + l5 )S 



(128) 



S R = -(l- l5 )S = S* L 



(129) 



so that Eq. (jl26J) generalizes to the superfield equation 




(130) 
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or, 



H^S L = -imSl (131) 
Now consider the supersymmetry transformation 

5S L = (r A a + S $ AP)S L = ~S L (132) 
This transformation commutes with the operator II^<9 M so that 

U%{S L + 5S L ) = m(S L + 5S L y. (133) 

If ipL is a left handed quark and B a (x) an antidiquark with the same mass 
as the quark, Eq. (jl32|) provides a relativistic form of the quark antidiquark 
symmetry which is in fact broken by the quark-diquark mass difference. The 
scalar supersymmetric potential is introduced through m — ► m + V s as 
before and Eq. ljlSOJ) remains supersymmetric. By means of this formalism, it 
is possible to reformulate the treatments given in the earlier sections in first 
order relativistic form. 

To write equations in the first order form, we consider V and $ given in 
terms of the boson fields by 

V = i-y^ - -<v (134) 

and 

$ = i 75 + (135) 
In the Majorana representation 

V* = -V, and $ = $* (136) 
We now define the left and right handed component fields by 

V L = ^—-^V (137) 



and 



1 + 

V r = —^V (138) 
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so that 



V R = -vi 

with 



-it .1-75 T/ I-75I T/ 



and making use of Eq. ()12U|) we have 



v L = in%- l -^ l -a, u v, v 

Noting that for any a M and b v we can write 



tt rr* u 1-75 1+75 h 1-75 
n^a^ = — - — 7^ — - — 7^/A = — ^ — 7 M 



and incorporating the definition of cr 



[IV 



I(n M n:-njip = l^^ 



in Eq.dHIJ) leads to 



Vt = iiijv„ + I(n„n; - n„ny 

Letting Ej„ = 11,11* - E^II*, Kq.(TPRj) reads 



We now have a first order equation: 



Similarly 

nj9 M y x = mV R 

Therefore 

nJnJ^^VL = mn;^ = m 2 v L 
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which after substitution of Eg . ([146)1 gives 

UV L = m 2 V L (149) 



For the $ part we can write 



with 



and 



nft^R = m$ L (150) 



$l = (151) 



*R = (152) 



Clearly, 



$ r = $t- (153) 
and above procedure can now be repeated for the $ fields. 

2 Skyrmion Model 

The Lagrangian density 



C = --F* tr(d„U l a"U) + tr([U^d„U, V^'OJJf) (154) 



consists of the original chiral model term (first term) plus the Skyrme (sec- 
ond) term. Here F n is the pion decay constant whose experimental value is 
186 MeV, and U is an 577(2) matrix given by 

U = exp{2ir ■ t^l} (155) 
Fix 

In above equations <f> is the pseudoscalar pion field (I = 1, S = 0), and e 
is an unspecified dimensionless parameter of the model. U is subject to the 
constraint 

UU ] = 1 (156) 
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and C is invariant under SU L (2) x SU R (2) under which 



U -> Lf/i? 1 " (157) 

where LL^ = 1, = 1 and r are the 2x2 Pauli matrices. F(r) is the 
shape function to be determined by the variation of the ground state energy 
of the soliton. In the large F w limit we can write 



so that 



which gives 



U ml + 2iT~ (158) 



d^d^U = 4F" 2 9 M • d»<f> (159) 



~ tr dprfdHJ = -\d^ ■ 0»4> (160) 

the kinetic term for a massless pion. 

A mass term for the pion would come from 

C m = tr (U + t/ f - 2) » -\ml4> 2 + 0(F n ) (161) 

which would break the symmetry of £ to the diagonal SU (2) (isospin) under 
which 

U -> MC/M t (L = R = M) (162) 

or 

r • -> Mr • 0M f (163) 

Going to the next order 

^ = 1 + 2<t.-^-2^ + ... (164) 



where 2 = • <j>, and 



7T 

we have 



U^ = l-2tr-^-2^ + ... (165) 



c^C/ = 2it ■ F-%<f> - 4F^ 2 cf> ■ d,<j> (166) 
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We can now write the right handed current 

J? = ^fU^d^ = ir-{<t>xd,4> + ^d,ct> + 0{F^) (167) 

where 

V^ = cf>xd^ + --- (168) 

is the isospin current, and 

A M = ^0^ + ... (169) 
is the axial vector current. Also 

■tf = Ud,U^ = -S- (d„U)tf = it ■ (V, - A,) (170) 

with additional terms coming from the Skyrme term. 
Under SU L {2) x SU R (2) we have 

U -> LUR\ U ] -> RU^l) (171) 

so that 

J M L - LJ^L\ J* - fljjflt (172) 

Both currents are conserved from the equation of motion (m T = 0). It is 
the charged components of that couple to the weak vector bosons 
in weak interactions. Both occur {Jb R ) in strong interactions governed by 
QCD. 

L is assumed to arise from the QCD theory of gluons and quarks when 
these colored fields are integrated out in the Feynman path integral to leave 
a skyrmion like effective theory in the color singlet sector. Note that using 

djfl = -U\d^U)U^ (173) 
C can be cast in the Sommerfield-Sugawara form 
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16 ; 32e 2 ^ J ; 

= ^ tr(J*J R ») + tr([J* J?} 2 ) (174) 
16 ; 32e 2 ^ v ; 

These currents obey the SUl(2) x SUr{2) current algebra of Gell-Mann, 
so that we have the correspondence with quark bilinears 

Ji~q l - Z ^i»q = v»-A, (175) 

and 



^7m9 = ^ + 4* ( 176 ) 
Another conserved current is the "topological" baryon current^ 

A /iua/3 

B" = ^ -pr tr{J?J*jK} (177) 

which is conserved identically without use of the equation of motion. The 
same current is obtained by using instead of J„. Integration over the 
baryonic current B° gives the baryonic charge. 

Lagrangian in Eq. (jl54|) should also contain the Wess-Zumino action 1 -' 
which can not be written in closed form in four dimensions but in five di- 
mensions 



wz : 



tr[a b ] (178) 

M 5 



240vr 

where the 1-form a is defined as a = Wd^U dx^ and M 5 is a five-dimensional 
manifold with boundary <9M 5 = M 4 , the Minkowski space. In the case of 
SU (2) the Wess-Zumino term vanishes. We can now proceed to show this. 

£ has extra symmetries that QCD does not have. For example, K + K~ — ► 
7r + 7r°7r~ violates parity P which is a symmetry of C. The Euler-Lagrange 
equation for C is 

dA F l U % U ) = ( 179 ) 
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We can try to add an extra term that violate P. Levi-Civita symbol t^ va i3 is 
a perfect choice. The Euler-Lagrange equation is generalized to 

d^Fl U^d^U) + Xe^ (U^d^U)(U^d u U)(U^d a U)(U^dpU) = (180) 

o 

We can now try to get the term in the Lagrangian that will yield the above 
equation of motion. In the equation above, the term that multiplies A is an 
obvious choice but it vanishes identically. To solve the problem, we write 
down the Wess-Zumino term which is a five- dimensional integral 

Cwz= I d 5 x tr[L^L u L a L^LA (181) 

JD5 

where L M = Wd^U, and D$ is a five-dim disc whose boundary M4 is space- 
time. Then C + Cwz would give the required equation of motion. 

Let Q be the generator of the gauge group. The U transforms as U — > 
U + ia(x)[Q, U\. Now introduce gauge field which transforms as 

-> - d^a(x) (182) 

L is not gauge invariant, so we have to add terms to it. The gauge invariant 
W-Z term is 

£\vz = £wz — J d 4 xA^J^ 

+ * J <f x ^ (d,A u )A a tr[Q 2 (d p U)lJi 

+ Q 2 U\dpU) + QUQU\d p U)U^\ (183) 

where 



J" = ^ tr[Q{d v Ulf){d a Ulf){d p UU^) 

+Q(U^d y U)(U ] d a U)(U ] d p U) (184) 

is the current coupled to gauge field A^. 
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At quark level, the baryon number operator is 




(185) 



Therefore, we have the baryonic current 

B " = ^ tr[{U^d v U){U^d a U){U^d p U)\ (186) 

Cyvz vanishes identically when there are only two flavors, i.e. U £ SU(2). 
This can be seen easily if we write 

U = l+ir-a (187) 

so that 

Cwz ~ J d 4 x ^ va<i aad^abdvacdadddpae tr[T a T b T c T d T e ] (188) 

In SU (2) there are only three independent generators r. However, the t ^ uaf5 
has four indices which require b, c, d and e be completely antisymmetric. 
This is impossible, so this forces £ = 0. The baryonic current, although it is 
obtained from Wess-Zumino term, is non-vanishing. 

The Skyrme Kink Solution; Rotating Skyrmion 
and its Quantization 

The equations of motion obtained for C have a static solution of the hedgehog 
form 

U = U (r) = exp{iF(r)^-} (189) 

r 

where r = x 2 + y 2 + z 2 , and r is aligned along isospin giving I + J invari- 
ance and 

U o (0) = -1, U (oo) = 1 (190) 
corresponds to the following boundary conditions on the function F(r): 

F(0) = tt, F(oo) = (191) 
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Here F(r) is the shape function to be determined by the variation of 
the ground state of the soliton. A realistic model would necessarily contain 
the dynamical quantities of spins and isospins. The time dependence of the 
soliton is obtained by rotating it with a time dependent 577(2) matrix A(t) 

U(r,t)=A(t)U Q (r)A\t) (192) 

where 

A(t)A\t) = 1 (193) 

Thus A(t) is a rotation matrix which makes the soliton rotate in time. 
Sustituting Eq. (|192j) into Eq. (jl54j) we obtain the Lagrangian 

L(t) = J d 3 x C(t, x) = -M + X tr(AA^) (194) 



where 



with 



• dA(t) . 2tt A , . 

at 3 e 6 F 1T 



A = / dp p 2 sin 2 F{\ + 4(F /2 + —^)) (196) 



sin 2 F , 

>o P 
Calculation gives A ~ 50.9 for the Skyrme solution. Now if 

A , • f a Q + ia 3 ia 1 + a 2 \ , . 

A = a + xa ■ t = . 197 

V ICL\ — Cl2 CLq — ICL3 I 



a 2 ) + a = 1), then the Hamiltonian is 



i=0 

with 

E a * = 1 ( 199 ) 



3 



i=0 
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constraints. This constrains a, on the sphere S 3 and (H — M) is the Laplacian 
on S 3 . Its eigenfunctions are spherical harmonics, polynomials in a ± ia 3 
and a± ± ia 2 of degree £. 

One obtains a tower of states with (spin=isospin) 

I = J = I (200) 

i = 1 and I = 3 are identified with iV and A respectively with the corre- 
sponding Ei eigenvalues of H giving their masses. 

The states for i > 3 are not found in nature. This is one of the weaknesses 
of the model. 

The tower can be chopped off by quantizing A(t) fermionically, so that 

{a l ,U j } = 5 lJ (201) 

In this case one has only i odd solutions. To have a non vanishing cubic 
polynomials in a« one has also to introduce a color degree of freedom for a« 
and choose color singlet cubic expressions for A. 

The wave functions are traceless symmetric polynomials in eij and that 
they have spin and isospin given by Eq. (|200j) can be immediately seen by 
considering the operators 

1 d d d 

J k = -«(-ao7 h a k - e k i m ai- — ) (202) 

2 oa k oa Q oa m 

and 

1 d d d 

h = -«(ao7; a-k^ e Mm ai- — ) (203) 

2 oak oa oa m 

The baryon charge density and the baryon currents are given by the 
components B° and B of the baryon current B 11 . They are 



and 



B l = — F> r jSk (205) 

7i z r 
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Note that the expressions for the spin and isospin are 



and 



with the baryon number 



K = -- tr(rA^A) (206) 



S = - tr(A^rA) (207) 



s= r^ w= mz£N! pes) 

Jo 7T 

Note that B = 1 when F(0) = tc and /(oo) = 0. 
We also note that integrating £ we obtain 

L = J d 3 r £ = -M + 2A5 2 = -M + 2\K 2 (209) 
where M is given by 



.sin 2 F ^ 1 sin 2 F sin 2 !-' 

and the moment of inertia A is 



M = 4* / rfr r 2 {^[F' 2 + 2^-] + ^ ^ + 2F' 2 }} (210) 



47T r°° F 2 2 S?77 2 F 

A = — / dr r 2 [^sin 2 F + ±stn 2 F(F' 2 + -^-)\ (211) 
3 y 2 e 2 r 1 

Let us now show these results explicitly. First we add a term ^f-B^B^ 
to C in Ea. dTolll . Since 

5^ = — — e^ Q/3 tr[U ] d v U U ] d a U U^dpU] (212) 

247T 

and writing U in component form, we have 

U = — A(cosF(r) + it ■ x sinF(r))A ] = — A(c + ir ■ x s)A ] (213) 
where we defined c = cosF(r), s = sinF(r) and D = det A. Now 

K a = -i tr(r a A ] A) = -i tr(Ar a A ] ) = -i tr(A^Ar a ) (214) 
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and 



S a = -i tr(A ] T a A) = -i tr{r a AA r ) = -i tr(AA ] r a ) (215) 
Therefore 

F 2 s 2 K 

where we defined 5j b = 5 a 5 — £ a :r fe . Integrating £ yields 
L = J £ d 3 r 

roo F 2 2 i 2 2 

2 V 2vr 2 r 2 ' 1 3 J 1 8 

= -M + 2A K • K = -M + 2\S ■ S 

= -M + \(K-K + S-S) (217) 

giving the desired expression. 

Now with the use of boundary conditions -F(O) = n and F(oo) = the 
shape function F(r) is solved numerically by minimizing M. The differential 
equation for F(p) is (where p is a dimensionless parameter defined as p = 

eF n r) 

p 2 (p 2 + 8sin 2 F)F" + 2p 3 F' + sin2F(4p 2 F' 2 - Asin 2 F - p 2 ) = (218) 
Behavior of F(p) is plotted against p in Fig.(l) below. 
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We further note that using Eq. (|154|) the equation of motion is written as 

VirfdpU) - -j— d»(tfd v U[rfd m uU, tfd v U\) = (219) 
and in the static limit we have 

diiU'diU) - J—d i (U^d J U[U^d l U, U^djU]) = (220) 

6 7T 

Substituting the soliton ansatz Eq. (|189|) into Eq. (j22(J|) we obtain 

r 4 F" + 2r 3 F' - r 2 sin2F + -^[sin2F(sin 2 F - r 2 F' 2 ) - 2r 2 sin 2 F F"\ = 

(221) 

which is identical to Eq. ()218p when we change the variable r to p = eF n r. 
therefore the hedgehog soliton is obtained from solving Eq. (j218|) is consistent 
with the above equation of motion. 

To find the mass spectrum of the solitons, we first have to write the 
system in the Hamiltonian form and then quantize the collective coordinates 
A. The matrix A has a canonical conjugate matrix A with matrix elements 

A a(i = dL/dA a(3 (222) 

where A a p is the matrix element of A. The quantization prescription is given 
by 

[A a p,An=i^ 

[A af3 ,A^] = [A a ^A^] = (223) 
The Hamiltonian is obtained from the Legendre transformation 

H = M + K 2 /2X = M + S 2 /2X (224) 

The baryon charge density and the baryon currents were given by the 
components B° and B of the baryon current above. From these the 
electric charge radius and the magnetic moments can be computed. 

Eigenvalues of the Hamiltonian are 

E = M+ — £(£ + 2) (225) 
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where i = 2 J. Therefore the nucleoli (N) and the A masses are given by 



MN = M+7 kl (226) 

Ma = M+ 2AT (22?) 
where M is obtained numerically from Eq. (j21()j) which comes out to be 

F 

M = 36.416 -j- (228) 
and similarly A is numerically evaluated from Eq. ()211|) to be 

A = y (-^r) x 50.42 (229) 

Above values are, for M = 865.15 MeV, and for A = 5.11 x 10~ 3 . Exper- 
imentally M N = 938.5 MeV, and M A = 1232 MeV. Best fit to these nucleon 
and delta masses yield e = 5.43 and F n = 129 MeV which is much smaller 
than the experimental value of F n = 186 MeV. Improvement of F„ value will 
be shown later in this paper. 

The isoscalar mean radius is 



< r 2 >e,i=o 



/*oo 2 /"OO 

/ 4nr 2 drBo(r)r 2 = / r 2 sin 2 F F' dr 

Jo 11 Jo 

2 f°° 4 47 

/ r 2 dr F' sm 2 F = (230) 

Jo « ' F 2 

gives 

< r 2 >J0 =O = 0.597 (231) 
The isoscalar magnetic mean radius is given by 



3 f n 4vrr 2 dr B (r) r 4 
< r 2 > M ,i=o= r ^ ~ 2 = 0.837fm 2 232 

5 <r 2 > E j=o 



so that < r 2 >m,i=o= 0.915 fm. 
Magnetic moments are given by 



2 ^ \ 

/ i p = 2M jV ( <r i2 ^ =0 +-) (233) 
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and 

/ i n = 2M JV ( <r i 2 2 ^ =0 -^) (234) 

Numerical calculations give (in terms of Bohr magneton) fi p = 1.871 and 
fi n = —1.34. The ratio |^-| = 1.4 is in excellent agreement with the experi- 
mental value 1.46. 

The axial coupling is given by 

9A = ~D (235) 

where 



f°° , r sin2F r . , 9 sin 2 F,, ,, 8sin 2 F,^ 
D= r 2 dr{—-[l + 4(F' 2 + — r -)] + F'(l + — — )} (236) 

Jo ' ' 

We find D = -17.2 and g A = 0.611. 
Strong couplings are given by 

2M N 

9-nNN = —^r- 9a (237) 

which is the Goldberger-Treiman relation. We get g^NN = 9.327 and g n NA = 
\g-nNN gives g w NA = 13.99. Experimental values are g^NN = 13.5 and g n NA = 
20.3. 

From the axial vector current and the vector current one can com- 
pute g w NN an d — . One finds 



9a = \ gT (238) 

and g n NN smaller than g^ N by 40%, as well as the magnetic moments ji p 
and fiN being smaller than the experimental values by 40%. We also fin 
|-^r| is good to 2%, and < r 2 > is off by 30%. Thus the agreement is only 

qualitative. Similar calculations performed by Nam and WorkmanEHl anc [ 
Hahm, et.al.^2 agree with these results. Magnetic moments are improved by 



adding vector mesons to the model'^', but F n resists all improved models so 
far. With many other contributing factors not yet taken into account, such 
as the effect of the pion mass and the interaction of vector mesons, this level 
of agreement is remarkable. 
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Possible Strategy for an Improved Skyrmion 
Model for Hadrons 

We have seen that the two-body potential models for the q — q and q — D 
bound state systems give a good description of hadrons when the poten- 
tials are extracted from QCD. In this case gluons are integrated out but the 
quark degrees of freedom remain. The model reflects the basic symmetries 
of QCD: its flavor independence, its approximate spin independence of the 
confining forces, the vector nature of gluons (Coulomb force at short dis- 
tances) and its chiral invariance broken by qq condensates. In addition we 
have q — q structure for mesons and qqq or q — D structure for baryons giv- 
ing the corresponding representations of £77(6) for low lying hadrons and its 
super symmetric extension SU (6/21). 

The skyrmion model, to be successful must reflect these symmetries. It 
does not. We now make the following observations toward the implementa- 
tion of such a programme: 

(1). Only I = J states are present. The approximate SU(A) symmetry 
of non strange hadrons has for the multiplet 7r, and (15 rep) and for 
the baryons N and A the 20 dimensional symmetric multiplet (submultiplet 
of 56 of SU(6)). This tells us that N and A are cubic functions of the quark 
wave function. 

In the skyrmion theory uj^ (I ^ J) is missing. One must decouple isospin 
from spin. For that one should transform the skyrmion rotation 



A(t)U (r)A\t), (AA^ = 1) 



(239) 



B(t)U (r)B\t), (BB^ = 1) 



(240) 



Uo(r) 

but also by deformation 
Uo(r) 

so that one has 



U(t,x) = M(t)U (r)N^(t) 



(241) 



with 

M(t) = B(t) A(t), N(t) = B\t) A(t). (242) 
Then one has 6 degrees of freedom instead of 4 and I and J can be decoupled. 
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(2). In the skyrmion theory wave functions of N are linear in A(t) and 
A is cubic in A(t). The SU{4) symmetry between N and A is lost. 

It can be recovered by having both N and A being cubic in M(t) and 
N(t). Then we can ask what do the I = 1 (linear in M and N) and the £ = 2 
(quadratic in M and N) represent? If one has the color degree of freedom, 
necessary for fermionic quantization, these states will be colored. Then it is 
natural to associate £ = 1 linear functions with the quark (cubic root of the 
Skyrme soliton) and £ = 2 with the diquark (|th root of the skyrmion). 

A simple modification of the skyrmion Lagrangian was proposed^^ with 
the aim of reflecting the cubic nature of the baryon. Uq(x) is replaced by 
Vq(x) = Uq(x) where 



so that 



U(t, x) = exp{2i t ■ (243) 



V(t, x) = exp{6i t ■ (244) 
The new Lagrangian CJ is adjusted so that for large it gives 

C! ~ - l - d„(t> ■ (245) 

It reads 

C " = tr(d,U^U n ) + — i-j tr([U^U n , U^ n d v U n } 2 ) (246) 

with n — 3. We have L = L\ for the usual skyrmion. Under SU{2) x 577(2) 
U n transforms as 

U n -> LU n B) (247) 

The hedgehog solution is obtained by having Uo(0) = —1 , Uq(oo) = 1 so 
that 

f/ 3 (0) = -1, f/ 3 (oo) = 1 (248) 

is unchanged. Again 

U (r) = exp {iF(r) (249) 
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with F(0) = Ti, F(oo) = 0. 

F(r) now satisfies a new linear equation derived from the new Lagrangian 
£3. The constants M and A have now new values as functions of F„ and e. 
We can now proceed to calculate the static properties of composite solitons. 
Rotating composite configuration is given by 

U n (r, t) = A{t)U%{r)A\t) (250) 
and its substitution into Eq. (j246|) gives the Lagrangian 



L n = -M N + 2X n s 2 (251) 
The mass of the composite static soliton is 



f 00 , 9r FiL ;9 sin 2 nF 1 1 sin 2 nF r sin 2 nF 
M n = 4vr dr r 2 {^[F' 2 + 2——] + — — — [— — - + 2F> 2 }} 



n 



n 2 r 2 2e 2 r 2 n 2 r 2 

(252) 

and the corresponding moment of inertia is 



4tt r° , 2 F 2 . 2 2 sin 2 nF l2 . sin 2 nF s 



A„ = — / rfr r^ S m 2 nF + — (2F' 2 + -^-)] (253) 

3 ./ 2rr e 2 n 2 n 2 r 2 



The normalized baryonic current 



= 2 2 e^ Q/3 tr[(t/-"^f/ n ) (U~ n d a U n ) (tr^CT 1 )] (254) 



has a time component 

R„ = - 

2vr 2 

and spatial component 



1 sin 2 nF 

B o = -S3 — — ( 255 ) 



B< = ^ ^ F< f, Sfc (256) 
7i z r 

The baryon number is equal to unity when F(0) = 71 and F(oo) = 0. 
Note that all odd powers of U obey the same "kink" boundary conditions 
U(0) = U 3 (0) = U 5 (0) = ■ • • = -1, and U(oo) = U 3 (oo) = U 3 (oo) = ■ ■ ■ = 1. 
This is not true for U 2 or other even powers of U. 
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The equation of motion for the shape function F(p) now takes the form 



p 2 (n 2 p 2 + 8sin 2 nF)F" + 2n 2 p A F' + -sin2nF(4n 2 p 2 F' 2 - 4sin 2 nF - n 2 p 2 ' 



n 





(257) 

Once again this is solved numerically and its behaviour is plotted for various 
odd values of n in figure 1 below. 



F(P) 
4 



\ \ 



11=1 



n = 3 



n = 5 



Fig.l Shape function F(p) vs. p 
Now the Hamiltonian is given by 

S 2 



H n = A a(S A aP -L n = M n + 



2\ r 



(258) 



In the usual manner, when they are fitted to the N and A masses one 
finds e = 1.72 (instead of 5.43 for n = 1), and F w = 185 MeV (instead of 129 
MeV for n = 1). 
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Isoscalar mean radius calculated in the same manner as before turns out 
< r 2 >]= = 0.74 fm, and the isoscalar magnetic mean radius is 1.185 fm. 
Similarly the magnetic moment calculation yields (where A in above [i v and 
fi n is now replaced by A n ) yields the values /i p = 2.0343 and ji n = —1.17. 
The ratio |^| = 1.97 (1.46 experimentally). 



The axial coupling is given by 



9a = ~ D n (259) 



where 



^ f°° 9 , r sin2nF r _ WW 9 sin 2 nF . 8sin 2 nF., 

Dn = l T dr{ ^~ [1 + 4(F +^V^) + F ( 1 + ^ 2 -)> ( 26 °) 

We find D 3 = —3.71 and hence g A = 1.31, which is close to the experi- 
mental value 1.23. The strong couplings for the n = 3 case turn out to be 
g^NN = 13.2 and g nNA = 19.8. 

Results of the calculations are compiled in the table below 



Quantity 


Old Prediction (n = 1) 


Prediction (n = 3) 


Experiment 


M N 


input 


input 


938.5 MeV 


M s 


input 


input 


1232 MeV 


F-n- 


129 


185 MeV 


186 MeV 


V< r 2 > E ,i=o 


0.597 


0.74 fm 


0.72 fm 


V< r 2 > M ,i=o 


0.915 


1.185 fm 


0.81 fm 


n P 


1.871 


2.034 


2.79 




-1.34 


-1.17 


-1.91 


9A 


0.611 


1.31 


1.23 


g-n-NN 


9.327 


13.2 


13.5 


9ttNA 


13.99 


19.8 


20.3 



We see that most everything is improved (some spectacularly). Except 
for the magnetic properties of the nucleon, the predictions derived from the 
model are in excellent agreement with the experimental values. Inclusion of 
the effect of vector mesons could change and improve the results significantly 
since the virtual photon can make a transition to the uj and p° mesons in the 
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vector dominance scenarios. Possible improvement of the predicted values of 
the magnetic properties will be discussed below. 

In the generalized model (as is also in n = 1 case) we find that the 
calculated nucleon charge densities oscillate within the nucleons, a peculiar 
future of nucleon substructure. It is well known, from the deep inelastic 
scattering experiments, that nucleons have point-like constituents. Since 
the quark degrees of freedom are integrated out in the Skyrme model, we 
expect the substructure of the nucleons not to be point-like. Therefore the 
oscillatory behavior of the charge densities may be due to interference effect 
of the constituents, which is similar to the quark interference effect found in 
the bag model^^. 

Excellent agreement with the experimental value of F n (=186 MeV) sug- 
gests that U 3 model could serve as a better description of the nucleons. The 
reason might be that composite model associated with the Lagrangian C n 
describes an effective hadronic theory based on QCD with the color group 
SU(n). The baryon consisting of n quarks are in a spatially symmetric 
ground state. In that case the current nB^ could be interpreted as a quark 
number current, so that the baryon has quark number 3 for n = 3. Under 
spin and isospin transformations, the skyrmion matrix 

U\r,t)=A(t)U^r) A\t) (261) 

has same transformation properties as the nucleon matrix of spin states n a 
and p a (a = 1, 2) 

N = ( ~ Hl Pl ) (262) 

with left and right multiplications being respectively associated with spin and 
isospin rotations. Now N and A which form the 20-dimensional symmetric 
representation of the spin-isospin group SU(4) are cubic functions of the 
elements of the underlying quark matrix 

(263) 

that can be taken to be proportional to A(t). When the color group is SU (3) 
and Q is quantized fermionically, it is possible to obtain the symmetric 20 
representation for the (N, A) system that is cubic in Q. Four of these twenty 
cubic functions, which are eigenstates of the Hamiltonian, are associated with 
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N. They have both spin and isospin 1/2 and can be put in a 2 x 2 matrix 
form QQQ, where QQ = Det Q with 

« - ( 2 -2 ) < 264 > 

being a singlet in spin and isospin. The remaining 16 cubic functions 
correspond to the wave functions of the A resonances. 

Since the quark has a baryon number 1/3 and the matrix Uq gives baryon 
number 1, it is natural to interpret the r dependence of the quark in the 
nucleon as being given by Uq{v). Then the quark matrix could be split into 
the form 

Q = A{t)U {r) (265) 

in a static approximation. The N and A wave functions would then be 
associated with the cubic functions of the elements of Uq as well as being 
cubic in A(t). The fermionic quantization is preferable for eliminating the 
higher excited states of the skyrmion represented by higher order polynomials 
in A(t) that are eigenstates of the Hamiltonian. A convincing derivation of 
n = 3 composite Skyrme Lagrangian from QCD theory of quarks and gluons 
with color group SU(3) is still lacking. For the moment one can regard it 
as a phenomenological model that gives reasonable quantitative predictions. 
We shall deal with this issue and the wave functions of the A resonances 
written in terms of the cubic functions in a forthcoming publication^^. For 
now we give a brief introduction to this subject. 

The baryon states can be written in terms of A, for example 

< A\A ++ , S 3 = ~ >~ A\ x (266) 

< A\p, S 3 = X - >~ A n (267) 

Note that we can write 

u |= a\ + »2 = ipx, d |= i(a,Q + ids) = ip2 (268) 

so that 

u i— — i(a — ia 3 ) — t/>2, d i— — (a x — ia 2 ) — — ipl (269) 
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giving 



with 



|^i| 2 + W 2 = l (271) 
Note that the left column of Q is 



while the right column is 



^ = -i(7 2 ^* (273) 
so that Q can be written in compact form as 



Now the A ++ state is 



and A + state is given by 



Q=(l> i) (274) 



A ++ = — ^ (275) 

7T 



A+ = — Vi(l - 3^ 2 ) = — ^i(2 - 3^i) (276) 

7T 7T 

In general we can construct the A resonances in this way. We can write 
down the QQQ (s — §) part in a matrix form 

qqq = ( -u\did 2 + u 2 d\ -u 1 u 2 d 1 + u\d 2 \ 
y u 2 did 2 — U\d\ u\u 2 d 2 — u\d\ J 

This term arises as part of a general quark matrix (in terms of it's and the 
gPs) as follows: Let us write down the following matrices 

QK a QK p Q + QKaQK a Q (278) 
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and 



QK a QKpQ + QKpQKaQ (279) 

each expression giving 3 x 4 — 2 = 10 triplet quark combinations as expected. 
In component form 

Ki = ^ K 2 = (280) 



and 



Rl = zi±* = (281) 



We can now write the full matrix Q as 



, -QKxQKxQ QK2QK2Q , 
1 -QK.QKO QK2QK2Q ' -~ ] 



where 



H,=-W,«=(^ ^fij (283) 

ni2 = QftQftQ=(_;^ (284) 

Si 21 = Q^Q^Q = ( (285) 

fi 22 = -QK,_QK,_Q = ( ""f ) (286) 

We can now embed s — § part into the center of the f2 matrix. We get 
the s — § part from 

Q\K2QK X + i^Qi^ - K 2 QK l - K 1 QK 2 ]Q = 

u 2 d\ - u^cfe -u 1 u 2 d 1 + u\d 2 . , . 
u 2 d\d2 — U\d\ uiu 2 d 2 — uldi 
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Similarly we can write down the s = | parts as 
Qi^QKi - K 1 QK 2 }Q = [ Ul d 2 - u 2 dt} 



and 



Q[K 1 QK 2 - K 2 QKi\Q = [u x d 2 - u 2 d x 



d x 
u 2 



Mi 

d 2 

Finally the full matrix containing all 20 terms become: 



(288) 



(289) 



/ d\ 



-U\d\ 



-u\di 



d\d 2 

—dido u 2 d\d 2 — Uid 2 



\ 



-d* 



Uid\d 2 + u 2 d{ 

2 



— U\ii 2 di + u{d 2 u{u 2 
U\U 2 d 2 — u\di —uiu\ 



u 2 d\ 



uid 



(290) 



2 "2 



The expression shown in Eq. 



then corresponds to the following term: 



—22 —23 
—32 —33 



(291) 



Note that the first column of H contains all the S3 = | A ++ , A + , A , 
and A - terms, the second column contains S3 = third column contains 



2 ■ 



S3 — and the last column contains S3 = expressions. 

We now give these terms explicitly For the A ++ states we have 

< A\ A++, S 3 = ~ >= ^ $ 

l 71 



<,si a + -. s. 4 = \>=-^ r 2 



< A\ A" 



5, 



< A| A++, S 3 
For the A + state we have 



-1 V6 
-3 _ V2 

Y >_ V 



^1 



^ 3 



3 v 7 ^ 2 



< A\ A + , S 3 = ->= — Vi ^2 

Z 71 



(292) 
(293) 
(294) 
(295) 

(296) 
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< A\ A+ S 3 = \ >= (1 - 3^2* ^2) (297) 

Z 7T 

< A| A+ S 3 = ^ >= — r 2 (1 - 3^! Vi) (298) 

Z 71 

< A\ A+, S 3 = ^ >= ^ ^ (29g) 

Z 71 

For the A states we have 

< A\ A , S 3 = § >= — V>i ^ (300) 

Z 7T 

< A| A , S 3 = i >= — ^2 (1 - ft (301) 

Z 7T 

< A| A , S 3 = ^ >= — Vi* (1 - M ^2) (302) 

Z 7T 

< A| A , S 3 = ^ >= — < ^2* (303) 

Z 7T 

Finally, for the A~ states we have 

< A\ A', S 3 = 3 ->=y^ ^ (304) 

Z 7T 

< A\ A-, S 3 = i >= V! ^ (305) 

Z 7T 

< A| A", S 3 = ^ >= ^ < ^2 (306) 

Z 7T 

<A| A" S 3 = ^>=-^^ 3 (307) 

Z 7T 



Action of the Hamiltonian operator given above on the A , A+, A~ and 
A ++ states will be dealt with in detail in a forthcoming publication. 
(3). Restoration of symmetries. 
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One can now have a supersymmetry algebra generated by the creation 
operators for I = 1 (fermionic) and those for i = 2 (bosonic). The I = 3 
cubic functions will be a representation of this superalgebra. 

(4). Construction of the enlarged Hamiltonian that will be approximately 
invariant under SU(4) and its supersymmetric extension 577(4/10) (subal- 
gebra of SU (6/21)) is the next step. 

Finding its eigenvalues and linear (quark), quadratic (diquark) and cubic 
(N and A) eigenstates is the following step. 

Final step will be recalculation of masses, charge radii, meson-nucleon 
coupling constants and magnetic moments, pion decay constant, meson- 
meson couplings, etc. 

Chiral Lagrangians with Vector Mesons 

In the large N c expansion baryons emerge as solitons within an effective 
theory of infinitely many weakly interacting mesons. A more realistic model 
should include not just the pseudoscalar meson field but also other heavier 
vector mesons. If we assume the 7r's, o/s and p's form an SU(4) multiplet, 
they should be treated in a unified manner. Accordingly, we now consider 
an 8 x 8 matrix configuration of these mesons 

W = exp{i lb T ■ 4> + -<J^ V + ^<vt ■ pT) (308) 

where 

d> — > — <b, uj — ► u, and p — > p (309) 

v F n ^ Ml ' H My K J 

will be substituted back into Eq. (j308|) later. Here g u and g p are the couplings, 
and M^ and M p are the masses of uj and p mesons. The 8x8 Pauli and 
Dirac matrices are formed by making the direct products 

r a -> r a ® 7 4 , and 7^ -> 7 2 ® 7 M 

also 

^ = ^[7^,7^] (310) 
First let us write down the simplest chiral model 
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C = tr[(d»W)W- 1 (d ll W)W- 1 ] + h.c. 

= -8 d»4> • + 4 (d^ aP + d„p af) ■ d»p a/3 ) + 0(4) (311) 

which is not sophisticated enough to contain much interesting physics. Since 
isospin space and spin space do not mix, there is no interaction among 
mesons. To connect the isospin space and the spin space while keeping the 
isospin symmetry intact we can insert r- and cr-matrices between the con- 
figuration W . To write down the more realistic Lagrangian we first evaluate 
the following terms: 

tr[y 9 (di i W)W~ 1 'y (d lt W)W- 1 + h.c] = 32 <9 M • d u (f> + 0(4) (312) 

and 

tr[ 1 ^W)W- l lv {d y W)W~ l ] + h.c. 

= 8 <9 M ■ d^cj> + 4 (dpUap + d,p ap ■ d»p af3 ) 

- 16 (d^ d x uj Xl/ + dpfT ■ d x p Xv ) 

- 48 dpfT ■ (0 x d v cj>) ~ 32 e^ K 4> ■ d» p{ d v u XK 

- 8 fff> ■ (d,p af3 x d vP n 

- 2 e, uSe e apXK p^ ■ (d»p XK x d v p 5€ ) + 0(4) (313) 

and 

tr [79 j IM (& l W)W~ 1 'f lu{d v W)W~ 1 } + h.c. 
= 16 • + % {dpUap d^ aP + d^p a(3 ■ d^p^) 

- 32 d^pT ■ (0 x d„<j>) - 64 e^p K 4> ■ d»p{ d u u XK 
+ 16 p a ? ■ (<9 M p a/3 x dupT) 

- 4 e^ e e a(3XK p a(3 ■ (d»p XK x d v p 5e ) + 0(4) (314) 
We can now combine Eqs. (j312l - I 314|) : 
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- tr [ le l^d^W" 1 ^ ^{d^W- 1 } + h.c. = 
= -8 d,cf> ■ 0"0 - 16 (d^ d x u Xu + d^pT ■ d x p Xv ) 

- 16 dppT ■ (0 x d v <t>) + 32 e^ K (f> ■ d» p{ d»uj XK 

- 24 pf* ■ (d»p af} x drpT) 

+ 2 e, uSt e a ^p al3 ■ (d»p XK x d v p Se ) + 0(4) (315) 
Now using Eq. (|3U9|) in W and Eq. (|315|) we have 

W = exp{^r ■ 4> + + ■ ( 316 ) 



and 



F 2 1 

-^■{tr[{d ,t W)W- 1 (d li W)W- 1 ] + -tr[7o(d"W)W- 1 7 (duW)W- 1 ] 
64 4 

+ tr[j^W)W- 1 -f u (d u W)W- 1 ] 

- tr [ le 7/ ,(cW)^-y lv {d v W)W~ l \ + h.c.} 

- jfc V • (0 x dv 4>) + F -g^ e^<j> ■ d»p{ Put" 

p u) p 



p 

+ <W e apXK p^ ■ (d»p XK x 9> fie ) + 0(4) (317) 

Since the masses of vector mesons are significant in comparison to the 
baryonic masses, we must include the meson masses in a more realistic La- 
grangian. Two mass terms comes from 



tr^W'fW - 4] + h.c. 

= 16 LO aP u aP + 16 p af3 ■ p afi - 16 e a(iXn 4> ■ p afi uj Xk + 0(4) (318) 
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which gives masses to the (spin triplets) vector mesons u and p while leaving 
the Ti massless (spin singlet). Using Eq. (j309j) we get 



F 2 M 2 

tr^W'fW - 4] + h.c. 



128 



F*M 2 g 2 t - F 2 M 2 g 2 



+ f :Z 9 m 9 ^ ^ ■ paP ^ + °w ^ 

p U) 



We now evaluate the term 



trlW- 1 n W n - 3] = -32 cp • </> + 16 p afi ■ p afi + 0(4) (320) 
and again using Eq. (|3U9|) we arrive at 

F 2 m 2 1 F 2 a 2 m 2 

triW" 1 n W r,-3] = -- ml 0-0+ ^ ' p a/? -p a/3 +0(4) (321) 

a term which shifts the masses of the isospin triplets it and p with respect 
to the isospin singlet u). We now combine all these into an approximate 
meson Lagrangian, expanded to third order, and containing pirn, ujpix and 
ppp interactions: 

^ d^pT ■ {cp x d v <p) + e^cp ■ d»p{ Put* 



M 2 ^ v ' M 2 M 2 

p U p 

o /?2 „3 | 

"J r IT JJp ntt /o . o U7,\ 1 

5 



8 ¥ 6- ^ • x - 2 m ' ' 



F 2 M 2 a 2 F 2 a 2 

e a px,4> ■ P Q/3 w A " + 0(4) (322) 
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4 M 2 M 2 

p a; 



P 

M 2 g p g w ^ M p Xfi 



Note that some of the total divergences have been subtracted. 
From the kinetic and mass terms of the Lagrangian we get 



Fl gl = 2Ml (323) 

and 

Fl gl = 2 M 2 p (324) 
which is the KSRF^l relation, and substituting 

M = M U) (325) 
we have mass relations for the mesons 

Ml-\ml = Ml (326) 

in perfect agreement with experiment. 
If we identify the u and p fields as 



dp = Ml uo u (327) 



and 



d, pT = M 2 p p v (328) 
up to second order, Lagrangian only of free fields (u mesons) 



1 i 



1 ^ 1 



d^-d x p Xu --p a ,-p^ (329) 
Take 

d^ v = Ml u", d^pT = M 2 p p v (330) 
So the equations of motion from C are 

d"u v - d v u» = lu" v (331) 
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□ = m\ 



(332) 
(333) 



Vector Mesons as Gauge Bosons 

Gauge potentials can be obtained from the decomposition of a larger algebra 
into a subalgebra and a coset. Gauge bosons can be combined with the 
fundamental mesons to form a single multiplet so that all mesons are treated 
equally in the larger symmetry group. 

We can write U in terms of the new variables E L and S r 

U = E L E R (334) 

Under the group 

SC/(2) L)global x SU(2) R , glohal x SU(2) local (335) 

El,r transform as 

E L -> 9(r) E L (336) 



E R - 9(r) E R fit (337) 

so that the Lagrangian 

C = W) (338) 

is invariant under the global symmetry A — > AU B\ with A G S , C/(2)x, )g iobai, 
£ G S£/(2) Rig i obal , and 0(r) G 5C/(2)i oca i. Since C/ transforms as C/ -> A C/ 
without any explicit 0(r) dependence, the symmetry is said to be "hidden". 
The associated gauge field 

V, = \ t ■ V, (339) 
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transforms as 



-> - 6(r) d M 9{r) + 9{r) 0\r) (340) 

We can now construct a Lagrangian that respects this gauge symmetry. Writ- 
ing the covariant derivative 

D lt = d /t -igV fl (341) 

where g is the gauge coupling, we can write two independent second order 
terms that are gauge invariant: 

£v = tr[{D,~ L ) S[ + (D,E R ) 4] 2 (342) 

£a = fr[p M H L ) S[ - 4] 2 (343) 

so that a combination of the form 

£ = C A + a C v (344) 

is equivalent to the non-linear sigma model Eq. (|338|) . Here a is an arbitrary 
parameter. Using unitary gauge 



we have 



~[=~ R =~ = m (345) 



Ca = tr(d^d^U) (346) 
16 

£v = tr{g V,- 1 - [(d,E) St + (0 M 5t) S]} 2 (347) 
We see that when we use the equation of motion for Cv 

g V; = -i tr{^ [(d,E) ^ + frtf) S]} (348) 

Cv vanishes and C becomes identical to Ca- Making the identification of 
vector fields with the p- meson fields p M , we add by hand the kinetic term 
of p, so that 



C = C A + a C v - \ Pilv ■ (f v (349) 
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Again, expanding to third term, using the weak field expansion we get the 
Lagrangian 




^ tr(d^ «W) --p F . fT 
g 2 Fl p„ • - | </ p„ • (p, x d» P(j> ) + ■■■ (350) 



a 



The mass of the p meson is given by 



(351) 



which is identical to earlier value when a = 2. Overall result of this approach 
agrees with experiment very well. 

Even though this description we just gave works very well with vector 
mesons, it would be more useful if we could find a way of including scalars 
and vectors in a unified manner. Toward achieving such a goal, we can 
now consider the genesis of covariant derivatives through a decomposition 
of a larger algebra involving spin or Dirac operators, into a subalgebra and 
a coset. The u-model obtained from the decomposition is valued in the 
subalgebra of the larger algebra, while the gauge fields are valued in the 
coset space. Once we identify the derivative of the coset elements as vector 
gauge potentials, the derivative of the larger algebra configuration can be 
written as a covariant derivative, i.e., a derivative of the subalgebra element 
plus a gauge potential. Let us rewrite W as 



W(<i>,u,p)=V(u,p) U(</>) 



(352) 



where 



U(<j>) = exp{— 75 r • (/)} 



(353) 



7T 



is a configuration of pion fields only, and 




(354) 
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U*V-%(VU) = tfV- 1 [Vd lx U + (d^V)U] 
U\d^ + V- 1 d,V)U = U^D^U 



(355) 



Therefore, the covariant derivative 



D 




(356) 



arises naturally with vector mesons being the gauge potentials V 1 d^ l V . 
Under this decomposition, the C of the non-linear a-model becomes 



£ ~ trlW-^d^W-^d^W)} -> tr[U\D ^JJ)^ (D^U)} (357) 



We see that the ordinary derivatives are replaced by the covariant derivatives, 
and the field configuration W is replaced by the field configuration U that is 
valued in the subalgebra. But there is one problem: we note that 



where F is a function of the vector meson fields, has a spin matrix a a p 
dependence but Wd^U has none. Therefore the trace of the matrix with 
scalar-vector meson interaction terms is zero. To generate non-zero values of 
such terms it is necessary to add a- and 7-matrices to C For simplicity, we 
now take 



V-%V ~ a a(3 d^F * 



(358) 




(359) 



D 




2M2 



(360) 



If we now define 



(361) 



and 



£<*/3 = V 1 (JapV 



(362) 
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such that they satisfy the same Dirac algebra as the original r- and a- 
matrices, namely 

{r a , I>} = {7 a , 7/3} = r] a p (363) 

and 

PajS,^] = V~ l [cr Q/3 ,cr^] V 

= 2i V" 1 (r] al jT,/3 U + rjp v E afl - Vav^pn ~ Vp^au) (364) 

Using these in Eq. ()322|) after some lengthy algebra we arrive at a weak field 
expansion of the Lagrangian 

£ = -\d^-d^-\ P , v -p" v 

- l - M 2 p p^pT- g p p„ • (0 x d»<p) + . . . (365) 

where we have made the identification d p p£ u = M 2 p p v m and obtained the 
KSRF relation. It is easy to see that this result agrees with hidden sym- 
metry approach we discussed. Therefore, a realistic connection between co- 
variant derivatives and the coset elements of the larger algebra has been 
demonstrated through this correspondence. 

Another example of obtaining covariant derivatives is from the consider- 
ation of a chiral Lagrangian that includes the meson nucleon interactions. 
The Lagrangian 

C = —ip 'jfj.d^ip ~ m ipWip + meson terms (366) 

is invariant under the chiral transformation. There is no mass m ipip for the 
fermions because such a term is not chiral invariant, hence the fermion field 
ip can not represent a realistic physical nucleon, which has to be massive. 
However, by performing a unitary transformation involving the meson field 
we can define a field £ 

£ = W i V = exp( Jr 7 5 r • + ^ a a(3 u^ + ^ a a/3 r ■ pl#) $ (367) 



66 



such that the new nucleoli mass term m becomes chiral invariant. Since 

W* = 74 74 = (368) 
we have £ = ^ Then, the transformed Lagrangian becomes 

£ = _(£ W-fy^iW-t - m & + • • • = -£7^ - m & + • • • (369) 
where 7^ = W'^'y^W'? and the covariant derivative = — with 

A, = = 75T • 9,0 + I 00, + I r ■ p M + • • • (370) 

We see that the Lagrangian in Eq. (j369|) gives a Dirac equation for the nu- 
cleons (7 -D M £ — rn £, which contains both the vector and the pseudoscalar 
meson-nucleon interactions through the covariant derivatives. 

These examples show that gauge potentials can be obtained from the 
decompositions of larger algebra into a subalgebra and a coset. Gauge bosons, 
which take on special roles in the model originally, can be combined with the 
fundamental mesons to form a single multiplet so that all mesons are treated 
equally in the larger symmetry group. 



Fermionic Quantization 

Earlier we quantized the skyrmion by rotating the soliton solution by the 
ST/ (2) time- dependent matrix A(t). The canonical conjugate matrix of A is 
given by 

A = dL/dA aP (371) 
and the Hamiltonian is obtained from 

H = A aP A aP - L (372) 

Although at the classical level A a p and A af3 are commuting variables, after 
quantization they become operators and are no longer commuting. We can 
either perform bosonic quantization 

[A aP ,A^]^i5% (373) 
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[A a p,A^\ = [A, A] =0 



(374) 



or fermionic quantization 



{A afi ,A afi } = i5iS; 



V 



(375) 



{A a/3 ,A^} = {A, A} = 



(376) 



that give the same algebra of the spin operators formed by the A's and 
A's. In the bosonic quantization, operators A a p commute among themselves, 
allowing particle states be described in terms of power functions in A of 
arbitrary orders. On the other hand, in the case of fermionic quantization, 
the A a p anticommute among themselves so that states with two or more 
identical A a p do not exist, and the powers in A of the state functions are 
limited. In the case when we take A to be fermionic, then the cubic functions 
for A ++ vanishes. 

To overcome this difficulty, we consider the same problem which occurred 
in the quark model when the problem arose of putting three spin up u quarks 
together to form s = 3/2 A ++ state. To conform with Pauli exclusion prin- 
ciple, three colors had to be introduced in the quark model to make the 
state antisymmetric in permuting any two quarks. In a similar sense, we 
now triple the number of variables in the Skyrme model by introducing two 
more matrices B and C in addition to the matrix A, and require the new 
Lagrangian to be invariant under the permutations of A, B and C. The 
color group SU(3) is replaced by the discrete subgroup Z 3 . In general for 
1/N expansion, the color group SU(N) is replaced by Zn, and there will be 
N different matrices. We propose a new Lagrangian that accommodates the 
new degrees of freedom 



We note that the time dependent parts of C remain unchanged, so we 
only have to concentrate on the time derivatives of 



C 



1 

3! 



{ 




tr{[U\ B d,U AB ,U\ B d u U AB f) 
B,C} (377) 
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Writing 



U AB {r,t) = A(t)U (r)B\t) 
U BC (r,t) = B(t)U (r)C\t) 

U CA (r,t) = C(t)Uo(r)A\t). (378) 



s A = - % - tr(rA^A) 



\ tr(TB*B) 



sc = ~\ tr(rC^C) (379) 

in the same way as before, after some tedious algebra, we find that the time- 
dependent parts of the non-linear sigma model terms are given by 



tr(d U AB d U AB ) = 2{s A - s B f - 8sinF cosF e abc s a A s b B r c 

+ 8sin 2 F(s A ■ s B — s A ■ r s B ■ r) (380) 



tr(d U BC d U BC ) = 2(s B - s c f - 8smF cosF e abc s a B s b c r c 

+ 8sin 2 F(s B ■ s c — s B ■ r s c ■ r) (381/ 



tr(d U CA d Ul A ) = 2{s c - s A f - 8sinF cosF e abc s a c s\ r c 

+ 8sin 2 F(s c ■ s A - s c ■ r s A - r) (382) 

The time-dependent part of the Skyrme term can also be similarly calculated. 
After some more lengthy algebra, we obtain the Lagrangian 



2F 2 ~ 

-M + --f V[(s A - s B ) 2 + (s B - s c ) 2 + (s c - s A ) 2 
6 lb 

2 

+ -X[(s A ■ s B ) + (s B ■ s c ) + (s c ■ s A )] (383) 
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where M and A are given by Eq. (|21Uj) and Eq. (j211j) . respectively, and V = 
f d 3 r is the size of the 3-space to be set equal to infinity after the Legendre 
transformation of the Hamiltonian. Defining the canonical conjugates 



E 



a/3 



dA aP 
dL 



dB, 



ai3 



dC a p 

and the corresponding spin operators 



S A = -~ tr(rA^A) 



= (384) 



2 

we obtain the Hamiltonian 



S c = ~ tr(rC ] C) (385) 



H = A al3 A a p + B a/3 B aj3 + C a/3 C aj3 - L 

= M + S V 6 US A - S B ) 2 + (S B - S c ) 2 + (S c - S A ) 2 } 
16 (QV - X) 2 

+{[(2F + X)S A + (2V - X)S B + (2V - X)S C ] x 

[{2V - \)S A + (2V + X)S B + {2V - \)S C ] 
+[(2V + \)S A + (2V - \)S B + (2V - \)Sc] x 

[{2V - \)S A + (2V - X)S B + {2V + X)S C ) 
+[(2V - X)S A + {2V + \)S B + (2V - X)S C ) x 

[(2V - \)S A + (2V -X)S B + (2V + Wo]} ~ _ (386) 
Finally, letting V — > oo, we arrive at the Hamiltonian 

H = M+^-(S A + S B + Sc) 2 (387) 
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whose eigenstates will be functions of the variables A, B and C . In analogy 
to the S77(4) quark model, we can describe the nucleons and the A's as 
SU (4) multiplets and write their states in terms of functions trilinear in A, 
B and C. For example, we can write the state function of spin 3/2 A ++ as 

<A,B,C\A ++ ,S 3 = ^)> ~ A 21 B 21 C 21 

+ antisymmetrization in A, B, C(388) 

and the spin up proton as 

<A,B,C\p,S 3 = ^)> ~ 2A 21 B 21 C 12 + 2A 12 B 21 C 21 + 2A 21 B 12 C 21 

—A\\B 22 C 2 \ — A 2 \B\\C 22 — A\\B 2 \C 22 
— A 22 B\\C 2 \ — A 2 \B 22 C\\ — A 22 B 2 \C\\ 
+ antisymmetrization in A, B, C (389) 

The problem of the infinite towers of particle states disappears because the 
Pauli exclusion principle does not allow any more states with I = S > 5/2. 

First Steps In the Implementation of SU (4) and 
577(4/10) Symmetric Skyrmion Program 

In a recent paper Cheung and Giirsey^^ showed that generators for SU (4) 
are constructed from dynamical skyrmion operators. In their paper they 
gave a fermionic realization of SU(2) x 577(2) algebra of the Skyrme model. 
They imbedded the SU{2) x S77(2) algebra into £7(4) by introducing new 
generators which are simple modifications of the SU (2) x SU(2) generators. 

The particle wave functions are written in powers of the collective coor- 
dinates. The infinite tower of states, which produces one of the problems 
for the Skyrme model can now be eliminated by restricting the values of 
these powers. We observe that if the collective coordinates are quantized 
fermionically, arbitrary functions of A (where tp ~ A p with 25 = P) can 
be written as polynomials in A of finite degrees, and therefore the height of 
the tower of particle states will have an upper bound. Another motivation 
for introducing fermionic quantization is that the predicted baryon masses 
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from the Skyrme model are too high if the experimental value of F n is taken 
as an input parameter. Modification of the model by addition of new terms 
will usually increase the particle mass further since usually the new terms 
give positive contributions to the Hamiltonian. Since in the renormalization 
procedure of a supersymmetric field theory the contributions coming from 
fermion loops are negative, such supersymmetric theory will have a better 
convergence. It is conceivable that the masses may be lowered if one can 
construct a supersymmetric Skyrme model, and fermionic quantization is a 
first step in this direction. 

Note that the spin and the isospin operators we introduced form the 6 
generators of the SU(2) x SU(2) algebra. To imbed the algebra into £7(4), 
we introduce the operators 

N = i tr(A ] A) (390) 

and 

E tj = 1 (KMV) (391) 

The form of the operators N and are easy to obtain as they differ from 
the operators S and K only by a factor of a Pauli matrix. One can easily 
verify that the 15 operators S, K and form the generators of S77(4). If 
we include the operator N, which commutes with all the 577(4) generators, 
we obtain a U (4) algebra. The SU (4) Casimir operator is given by 

C = S 2 + K 2 + EijE ij (392) 

The fermionic quantization prescription allows us to expand the algebra 
into superalgebra. The supersymmetric extension of the U (N) algebra and 
the operators forming the SU (2) x SU (2) subalgebra can now be obtained. 
Let A a be an TV-vector, with A a its canonical cojugate (a = 1,...,N). 
Together with their complex conjugates they obey the following commutation 
relations: 

{A a , Ap} = {A* a , A*A = {A a , Ap} = {A* a , A*A = (393) 
{A a ,A;} = {A a ,A}} = (394) 
{A a ,A* f ,} = {A* a ,A ft } = (395) 
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{A a , A/3} = Sa/3, {A* a , A*p} = 5 a/3 

Defining the operators J a p as 



(396) 



J af , = A a A f) + AZA* f) (397) 
They generate the U (N) algebra 

[J a /3, Jnu] = SpuJav ~ ^avJ^fi (398) 

For the case of N = 4, we keep the track of the imbedding by noting that 
the operators 

So = ^(Ju + J22 - J33 ~ Ju) (399) 
S + = -L(J n + J 24 ) (400) 
S_ = -L(J 31 + J 42 ) (401) 

h = \{Ju ~ J22 + J33 ~ Ju) (402) 
J + = -L(J 12 + J 34 ) (403) 

/- = -J=(J 2 i + J 4 3) (404) 

form an SU (2) x SU (2) subalgebra of the Skyrme model. 

So far we have only the even algebraic operators J a p that are products of 
two fermionic operators A and A. To make a supersymmetric extension of the 
algebra, we need some odd algebraic operators that are made of products of 
odd number of fermionic operators. The even and odd operators are expected 
to form a superalgebra. 

To obtain the superalgebra we now introduce the (odd) supercharge op- 
erators 



q 7 = xa, + ^a;a;a; (405) 
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and 



Q x = \*A* X + -A V A\A 



(406) 



where A is an arbitrary complex parameter. Using the commutation relations 
we see that J a p, Q 7 and Q 7 with (a, (3, 7 = 1, . . . , N) form the superalgebra 



{Q 7 ,Qa} = {Q 7 ,Q a } = 



(407) 



{Q 7 , Q\\ — J-yX ~ 8-y\Jv 



(408) 



[JaPi Q-y] — ~8a-yQ/3 

The quadratic invariant of the superalgebra is 



Co 



J 2 — J^fJuu ~~ JapJpa + [Qjj Q 7 ] 

which can be obtained by taking the supertrace of Q 2 , where 



(409) 



(410) 



n 



J 



IN 



Qi 



\ 



Qn 
T N J 



(411) 



Jni ■ ■ ■ Jnn 

\ ~Qi ■ ■ ■ ~~Qn 

is a supertraceless matrix configuration with even and odd operator matrix 
elements. Other invariants of the superalgebra are C3 = Str(Q 3 ), . . . , Cn = 
Str(n N ). 

Among the irreducible representations of SU(4) C U(4), there is a 15- 
dimensional boson multiplet (tt, p, u) and a 20-dimensional fermion multiplet 
(N, A). With the supersymmetric extension of U(4) it is possible to extend 
Lagrangians which contain only the tt, p (and uo) to the ones that contain 
mesons, nucleons and A's. The new Lagrangian must necessarily contain the 
Dirac and Rarita-Schwinger Lagrangians for iV's and A's as well as the usual 
meson Lagrangian. The Hamiltonian to be obtained then should resemble 
the casimir invariant C2 shown above, and the supersymmetry breaking will 
bring the Hamiltonian to the form Eq. (j392J) . and the subsequent breaking of 
£7(4) into SU(2) x SU(2) will bring the Hamiltonian to the form of En. 1(223)1 . 

In order to avoid the problem of predicting too high soliton masses when 
taking physical value of as an input parameter, many people subtract 
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from the Hamiltonian a zero point energy without explaining its origin. We 
suggest that the subtraction of the zero point energy without explaining its 
origin. We suggest that the subtraction of the zero point energy may come 
from the remnant effect of supersymmetry after the supersymmetric theory 
has been broken down to a non-supersymmetric theory. It is known that in a 
field theory where supersymmetry is exact, the ground state energy is zero. 
The positive contribution from bosons is exactly canceled by the negative 
contributions from the fermions. In a theory with approximate hadronic 
supersymmetry, although we do not expect exact cancellation of bosonic 
and fermionic parts of the soliton energy, it is probable that there is a good 
balance of both sides so that the baryonic masses can be lowered significantly 
to give a reasonable value of the pion decay constant. 

Conclusion and remarks 

Quark model with potentials derived from QCD, including the quark-diquark 
model for excited hadrons gave mass formulae in very good agreement with 
experiment and went a long way in explaining the approximate symme- 
tries and supersymmetries of the hadronic spectrum, including the symmetry 
breaking mechanism. Near success of the minimal SU (2) Skyrme model was 
modified to a generalized composite Skyrme model in which baryons are con- 
sidered as composite objects made up of many SU (2) configurations. In the 
special case of three SU (2) configuration, the static properties of the soliton 
had excellent agreement with the experimental data. 

According to the large N c expansion, a more realistic model should include 
the vector mesons and the pions. Inclusion of low lying vector mesons u and 
p in a modified model, where the Lagrangian of the model was constructed by 
inserting 7- and a- matrices that allowed us to get the usual meson physics, 
including the meson mass relations in perfect agreement with experiment. 
All the mesons were treated on equal footing so that they belong to a single 
meson multiplet. 

The simplified model we built, in which a soliton ansatz to the Lagrangian 
£ has not yet been found, gives us a good indication of the salient features 
of any soliton solutions to the full model. 

Magnetic properties are improved by adding vector mesons to the model. 
With the static approximation, the magnetic part of the configuration van- 
ishes and the model is then becomes reducible to an SU(4) Skyrme model 
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that produces a variety of mass spectra for the solitons. 

Fermionic quantization for the collective coordinates of the skyrmion is 
introduced. It has the advantage that through the Pauli exclusion principle it 
it chops off the unphysical infinite tower of (S — I > |) particle states. With 
this approach, nucleons and A's can all be expressed as trilinear functions of 
the collective coordinates as we have shown in this paper, and the fermionic 
quantization approach allows a simple incorporation of supersymmetry into 
the model. Such a supersymmetric extension makes it possible to lower 
the soliton mass since fermions contribute negatively to the Hamiltonian. 
Imbedding of the SU(2) x 577(2) symmetry into U(4) and U{A/n) (n = 1,2 
in our case) is the first step towards a supersymmetric extension. U(N) 
algebra is extended to a superalgebra for arbitrary iV by the introduction of 
supercharges, and the Casimir invariants are easily obtained by taking the 
supertrace of the powers of the configuration Q. 

The fermionic degrees of freedom must be introduced in such a way that 
an approximate symmetry under S77(4/10) holds. It is to be seen if after the 
suggested modifications one will obtain a quantitatively viable skyrmionic 
description of hadronic physics at least as good as the 2-body potential 
model 1 ^' 1 ^ inspired by QCD and reflecting its fundamental symmetries. 
It is highly desirable to find an explicit relationship between the relativis- 
ts quark model we built and the Skyrme model. These two models should 
have common hadronic symmetries. The approximate S77(6) symmetry and 
SU (6/21) supersymmetry have been established for the quark model, but not 
yet for the Skyrme model. We can build a supersymmetric Skyrme model 
based on the superalgebra we developed. Some preliminary work has been 
done in that direction, which will be the subject of another publication. It 
will be very interesting to compare the two different supersymmetry realiza- 
tions of the quark model and the Skyrme model. 

Strangeness should also be incorporated into the model. If we include 
strangeness in our model while treating all mesons equally, we would expect 
the mesons to form an SU(6) multiplet. With such modifications we would 
expect that the predictions of the static properties of the solitons can be 
brought into much better agreement with experiment. 

Also, it should in principle be possible to construct Fermi fields from 
Bose fields in 3 + 1 dimensions in the light of these recent developments. The 
picture that one would like to see emerge is a duality between a field theory 
written in terms of quark and gluon fields that would accurately reflect the 
hilbert space structure at very small distances on one hand, and a dual theory 
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expressed in terms of non-linear meson couplings whose soliton solutions 
would correspond to the structure of the Hilbert space as seen from a larger 
distance perspective on the other. This would require two boson-fermion 
connections. First the quarks and gluons must give rise to the non-linear 
meson fields as composites and these latter must be reinterpreted in terms 
of baryons. It is not clear how fundamental this phenomenologically hoped 
for connection can be. 

Other applications such as when one considers the pion fluctuations around 
a soliton one can study the meson-nucleon scatterings. In the many-skyrmion 
system, one can study dibaryon states, nucleon-nucleon interactions and 
properties of nuclei. Skyrme model is a viable model of baryons. Recent 
developments of the model give strong support to the belief that baryons can 
be regarded as topological solitons of a theory of interacting meson fields, 
and the low energy physics should be described non-perturbatively by an 
effective Lagrangian based on the underlying symmetries of QCD. 
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